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SECTION  I 


INTRODUCTION 

The  principal  purpose  of  this  investigation  is  to  provide  a  basis  for 
assessing  the  capability  of  the  computer  program  CLAPP  [l]  to  predict 
buckling  loads  for  fiber-reinforced,  circular  cylindrical  panels  under 
prescribed  uniform  axial  end-displacements. 

Two  testing  programs  were  undertaken  to  accomplish  this  objective. 

They  are  designated  as  Program  A  and  Program  B  in  this  report. 

Program  A.  Experimental  buckling  loads  were  determined  for  20  specimens 
having  identical  fiber  patterns  of  [0/90]^  and  20  specimens  having  iden¬ 
tical  fiber  patterns  of  [0/+45/90] ^ .  The  specimens  of  a  particular  fiber 
pattern  were  characterized  geometrically  by  five  different  aspect  ratios 
and  by  two  different  sets  of  boundary  conditions.  The  aspect  ratios  were 
a/b  =  1/2,  3/4,  1,  4/3,  and  2,  where  a  and  b  are  the  dimensions  of  the 
projection  of  a  panel  into  its  base  plane  as  shown  in  Figure  1.  Boundary 
conditions  along  the  straight  edges  of  a  specimen  corresponded  to  either 
an  unsupported  edge  or  a  simply-supported  edge.  For  the  simply-supported 
edge  a  distinction  between  unconstrained  and  constrained  circumferential 
displacements  is  also  made.  Boundary  conditions  along  the  curved  edges 
corresponded  to  a  clamped  edge  for  all  specimens. 

Program  B.  Experimental  buckling  loads  were  determined  for  6  specimens 
with  the  same  fiber  pattern  ( [0/90 1 2  s )  *  t*1e  same  aspect  ratio  (a/b  =  1.247, 
a  =  16  in.  and  b  =  12.83  in.),  and  the  same  boundary  conditions  (unsupported 
straight  edges  and  clamped  curved  edges).  Similarly,  experimental  buckling 
loads  were  obtained  for  5  specimens  with  the  fiber  pattern  [0/+45/90]s, 
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aspect  ratio  a/b  =  4/3  (a  =  16  in.  and  b  =  12  in.),  and  boundary  conditions 
consisting  of  simply-supported  straight  edges  and  clamped  curved  edges. 

A  second  purpose  of  this  investigation  is  to  modify  the  computer  pro¬ 
gram  CLAPP  so  as  to  minimize  the  effort  required  to  input  necessary  data. 
This  phase  of  this  investigation  is  contained  as  a  USER'S  MANUAL  in 
Appendix  D. 

Finally,  a  third  purpose  is  to  modify  CLAPP  so  that  it  can  be  used  to 
predict  the  buckling  behavior  of  f i ber- re i nforced ,  circular  cylindrical 
panels  and  flat  plates  that  are  augmented  by  longitudinal  stiffeners.  This 
is  accomplished  for  isotropic  stiffeners  with  a  variety  of  cross  sections 
and  for  a  quas i - i sotropic  stiffener  with  a  hat  shaped  cross  section. 
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SECTION  I  I 


TESTING  PROCEDURE 


(ll-l).  TEST  SPECIMENS.  Each  of  the  circular  cylindrical  specimens  of 
Program  A  and  Program  B  was  laminated  from  graphite-epoxy,  and  each  speci¬ 
men  was  cured  in  a  mold  for  which  the  external  radius  was  12  inches. ^ 

The  external  radius  of  a  specimen  tended  to  be  slightly  less  than  12  inches 
upon  removal  from  the  mold.  The  thickness  of  a  test  specimen  was  taken  as 
an  average  of  the  thicknesses  at  32  locations  along  a  perimeter  located 
uniformly  1-1/2  inches  from  the  exterior  perimeter  of  the  specimen.  The 
average  thickness  for  each  specimen  of  Program  A  is  shown  in  column  2  of 
Table  1.  The  average  thickness  for  each  specimen  of  Program  B  is  shown 
in  column  2  of  Table  2.  Figure  1 
indicates  pertinent  geometrical 
parameters  of  a  typical  specimen. 

The  radius,  R=  12.0  in.,  is  the 
radius  of  a  perfectly  circular 
panel  installed  in  the  testing 
fixture,  and  the  thickness  shown 
is  the  average  thickness  of  a 
typical  panel . 

Laminate  stiffness  properties 
were  calculated  by  lamination 
theory  [2]  using  the  following 
lamina  properties: 


T 


0.038  in. 


R=12.0  in. 


Figure  1.  GEOMETRY  OF  CURVED 
PANEL  SPECIMEN 


(l)  All  specimens  were  fabricated  and  cut  by  the  Air  Force  Flight  Dynamics 
Laboratory,  Air  Force  Wright  Aeronautical  Laboratories,  Air  Force  Systems 
Command,  Wr ight-Patterson  Air  Force  Base,  Ohio. 


=  752  ksi  ,  i 


=  0.335,  and  v 


E,,  =  20,524  ksi ,  E22  =  1,333  ksi  ,  G)2 
0.022.  and  E22  are  Young's  moduli 
dicular  to  the  fiber  direction,  G^2  is 
with  these  two  directions,  and  j  and 
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of  elasticity  parallel  and  perpen- 
th.  shearing  modulus  associated 
v.,  are  Poisson  ratios. 


Program  A.  The  specimens  of  Program  A  were  distinguished  by  aspect  ratio 
(a/b) ,  laminate  fiber  pattern,  and  boundary  conditions. 

Specimens  having  five  different  aspect  ratios  (a/b  =  1/2,  3/4,  1,  4/3, 
and  2)  were  tested.  It  is  convenient  to  perceive  of  specimens  that  are 
associated  with  a  specific  aspect  ratio  to  form  a  group.  Thus,  five 
separate  groups  of  specimens  are  identifiable. 

Each  group  of  specimens  consisted  of  two  sub-groups:  one  sub-group 
being  associated  with  simply-supported  straight  edges,  and  one  sub-group 
being  associated  with  unsupported  straight  edges.  Boundary  conditions 
along  the  curved  edges  were  clamped  for  all  specimens  in  each  group. 

Two  laminate  f  iber-patterns  ( [  0/90  3  2  s  anc*  1 0/+45/90  ]  )  were  considered 
for  each  sub-group.  Moreover,  each  member  of  each  group  was  duplicated  to 
provide  a  means  to  assess  repeat ib i 1 i ty  of  experimental  results. 

From  the  foregoing  discussion  it  will  be  observed  that  each  group  con¬ 
tained  eight  specimens:  four  with  simply-supported  straight  edges  and 
four  with  unsupported  straight  edges.  Of  the  four  specimens  in  a  sub¬ 
group  (either  simply-supported  or  unsupported  straight  edges)  two  were 
characterized  by  a  f iber-pattern  of  [0/90]^  and  two  were  characterized 
by  a  fiber-pattern  of  [0/+45/90]s.  The  complete  experimental  effort 
amounted  to  forty  specimens.  Table  1  has  been  organized  to  ref’ect  the 
foregoing  classification  of  specimens. 
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The  physical  dimensions  of  the  base-planes  (projected  area  of  a  specimen) 
of  the  specimens  associated  with  the  five  aspect  ratios  were  8  x  16,  12  x  16, 


16  x  1 6 ,  16  x  12,  and  16  x  3  (all  dimensions  are  in  inches).  These  dimen¬ 
sions  correspond  to  the  inside  dimensions  of  the  specimens  after  installation 
in  the  testing  fixture. 

Program  B.  As  was  stated  in  the  introduction,  the  specimens  of  Program  B 
can  be  categorized  into  two  groups:  one  group  with  fiber  pattern  [ 0/90 1 2 s ’ 
aspect  ratio  16/12.83  s  1 .247,  and  unsupported  straight  edges;  and  one 
group  with  fiber  pattern  [0/+^5/S0]^  aspect  ratio  16/12  =  4/3,  and  simply- 
suported  straight  edges.  The  results  of  buckling  tests  on  these  specimen 
are  tabulated  in  Table  2. 

(IJ-2).  TEST  FIXTURE.  The  testing  fixture  used  in  both  testing  programs 
was  a  modification  of  a  test  fixture  used  by  Wilkins  [3].  Modifications 
of  the  Wilkins'  testing  fixture  were  necessary  to  accommodate  specimens  of 
different  aspect  ratios.  Figure  2 
shows  the  head-plate  and  the  base¬ 
plate  (with  their  auxilliary 
pressure  blocks  and  clamping 
plates)  through  which  an  axial 
compressive  load  was  applied  to 
a  specimen.  All  surfaces  of  the 
head-plate,  and  of  the  base¬ 
plate,  were  carefully  machined 
so  that  they  were  within  0.001  in. 
of  being  paral lei . 


FIGURE  2.  Head-plate  and  base-plate  with 
auxiliary  pressure  blocks  and 
clamping  plates. 


TABLE  I.  EaPERIMENTAL  AND  NUMERICAL  BUCKLING  RESULTS  FOR 
CIRCULAR  CYLINDRICAL  PANELS  FOR  TESTING  PROGRAM 


(11-3)  .  TESTING  MACHINE.  An  axial  compressive  load  was  applied  to  j  speci" 
men  through  a  120,000  lb  Tinius-Olsen  hydraulic  testing  machine.  Most  experi¬ 
mental  buckling  loads  were  of  such  a  magnitude  that  the  intermediate  range 
(12,000  lb  range)  of  the  testing  machine  could  be  used  effectively.  The 
finest  division  for  this  range  is  50  lb.  The  low  range  (3000  lb  range)  was 
used  to  test  several  specimens  with  unsupported  straight  edges.  The  finest 
division  for  this  ranye  is  5  lb. 

The  surfaces  of  the  plat  ten  and  the  cross-head  of  the  testing  machine 
were  dressed  so  that  they  were  nearly  plane. 

The  head-plate  and  the  base-plate  of  the  testing  fixture  are  shown 
installed  in  the  Tini us-01 sen  testing  machine  in  Figure  3- 


Early  in  testing 


Program  A  i t  became 
apparent  that  the 
cross-head  would 
tilt  as  the  react i ve 
force  of  the  speci¬ 
men  on  the  cross-head 
t ended  to  lift  it  off 
the  threads  of  the 
vertical  columns.  To 
eliminate  the  tilting 
act  ion  two  I  a  rge 


aluminum  nuts  were 
machined  and  placed 
on  the  vertical  screws 


FIGURE  3-  Testing  lixture  installed  in  the  Tinius- 
Olsen  hydraulic  testing  machine. 
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beneath  the  cross-head.  At  an  appropriate  point  during  the  installation 
of  a  specimen  these  nuts  were  tightened  against  the  lower  surfaces  of  the 
cross-head  causing  it  to  lock  in  place.  These  aluminum  locking  nuts  are 
shown  in  Figure  3- 


(11-4).  IMPERFECTION  DEVICE.  Figure  4  shows  the  mechanical  device  that 
was  used  to  measure  deviations  of  a  specimen  from  a  perfect  cylindrical 
form;  that  is,  to  measure  the  initial  geometric  imperfections  of  a  speci¬ 
men  . 

Components  of  this  device  were  constructed  so  that  the  vertical  sides 
of  the  circular  groove  in  the  platform  of  the  device  were  concentric  with 


the  vertical  sides  of  the  circular  groove  in  the  base-plate  of  the  testing 


fixture.  The  imperfection  device 
was  secured  to  the  base-plate  of 
the  testing  fixture  so  that  its 
slotted  vertical  platform- 
supports  were  perpendicular 
to  the  base-plate.  This  en¬ 
sured  concentricity  of  the 
platform  groove  with  the  groove 
in  the  base-plate  for  any 
platform  leve I . 

To  position  the  imperfec¬ 
tion  device  relative  to  the 
base-plate  of  the  testing 
fixture,  the  dial  indicator 


mechanism  shown  in  Figure  4  cirnoc  <•  , 

FIGURE  4.  Imperfection  measuring  device 

attached  to  t he ' base-p I  at e  of 

the  testing  fixture. 
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was  designed  so  that  it  could  be  moved  smoothiy  and  snugly  in  the  circular 
groove  of  the  platform  of  the  imperfection  measuring  device.  With  the  tip 
of  the  dial  indicator  extension  resting  against  the  vertical  side  of  the 
circular  groove  in  the  base-plate,  the  mechanism  was  moved  along  the  plat¬ 
form  groove.  By  trial  the  imperfection  measuring  device  was  adjusted  rela¬ 
tive  to  the  base-plate  so  that  the  pointer  of  the  dial  indicator  was  undis¬ 
turbed  for  a  complete  cycle  along  the  platform  groove.  Once  the  proper 
position  was  located  the  imperfection  measuring  device  was  secured  to  the 
base-plate  of  the  testing  fixture  by  machine  bolts.  Positioning  pins  were 
installed  so  that  the  correct  position  could  be  duplicated.  This  procedure 
resulted  in  a  variation  of  less  than  0.001  in.  for  a  complete  traverse  of 
the  platform  groove. 

Holes  for  positioning  pins  were  drilled  at  one- inch  intervals  along  the 
slotted  vertical  platform  supports.  This  ensured  that  the  position  of  the 
imperfection  measuring 
device  relative  to  the 
base-plate  of  the  testing 
fixture  could  be  duplicated 
at  appropriate  levels. 

Once  the  platform  was  lo¬ 
cated  at  a  specified  level 
by  the  positioning  pins, 
i t  was  locked  i n  that 
position  by  cap  screws. 

Figure  5  shows  the 

imperfection  measuring  FIGURE  5-  Imperfection  measuring  device  in 

position  to  measure  initial  qeometric 
imperfections  of  a  specimen. 


device  installed  on  the  testing  fixture.  Deviation  of  a  specimen  from  per¬ 
fect  circularity  were  measured  by  movinq  the  dial  indicator  mechanism  along 
the  platform  groove. 


(I  I  -  5 ) -  POSITIONING  OF  HEAD-PLATE  AND  BASE-PLATE.  Two  requirements  were 
essential  for  proper  alignment  of  the  head-plate  relative  to  the  base¬ 
plate  of  the  testing  fixture.  First,  the  circular  arc  associated  with  an 
edge  of  the  circular  groove  in  the  head-plate  and  the  correspond i nq  cir¬ 
cular  arc  of  the  circular  groove  in  the  base-plate  must  lie  in  parallel 
planes.  This  requirement  is  realize"'  w‘- ->n  the  head-plate  (when  attached 
to  the  cross-head  of  the  testirg  r.,  •  and  the  base-plate  (when  attached 

to  the  plat  ten  of  the  testing  r  c  are  parallel.  Secondly,  these  two 

arcs  must  lie  in  a  conmon  c'rc  lindrical  surface  that  is  perpendi¬ 

cular  to  the  base-plate  (and,  hence,  pe rpend i cu 1 ar  to  the  head-plate). 


Paral lei  ism  of  the  head- 
plate  and  the  base-plate  was 
assessed  using  the  device 
shown  in  Figure  6.  This 
device  consists  of  a  rigid 
base  that  was  machined  to 


fit  snugly  in  the  circular 
groove  of  the  base-plate, 
and  that  could  be  moved 
smoothly  along  the  groove. 
A  dial  indicator  was 


a  t  t  ached  to  a  stiff  t  i'o  I 
rod  which  was  affixed  to 


FIGURE  f> .  Device  used  to  establish  the  relative 
pa ra 1  I e 1  i sm  of  the  head-plate  and  the 
base-plate  of  t  lie  test  inn  lixfnre. 


I? 


the  moveable  base.  With  the  plunger  of  the  dial  indicator  resting  against 
the  surface  of  the  circular  groove  in  the  head-plate,  the  device  was  moved 
along  the  circular  groove  in  the  base-plate.  Using  the  centerline  of  the 
circular  groove  as  a  reference,  it  was  observed  that  the  variation  in  the 
distance  between  the  surfaces  of  the  grooves  in  the  head-plate  and  the 
base-plate  did  not  exceed  0.0035  in.  on  either  side  of  the  centerline  for 
the  maximum  arc  of  1 8 . 5  in.  (ci rcumferential  arc  length  of  the  widest  speci¬ 
men).  This  variation  was  accordingly  smaller  for  narrower  specimens.  Once 
the  head-plate  and  the  base-plate  had  been  assessed  to  be  as  parallel  as 
reasonable  efforts  would  allow,  they  were  securely  clamped  to  the  cross¬ 
head  and  to  the  platten  of  the  Tinius-Olsen  testing  machine  by  specially 
designed  clamps.  These  clamps  can  be  seen  in  Figure  5- 

The  imperfection  device  was  used  to  bring  corresponding  circular  arcs 
in  the  head-plate  and  in  the  base-plate  into  concentricity.  To  accomplish 
this  alignment  the  platform  of  the  imperfection  measuring  device  was  posi¬ 
tioned  at  an  appropriate  level  and  the  tip  of  the  dial  indicator  extension 
shown  in  Figure  was  allowed  to  rest  against  the  vertical  side  of  the 
circular  groove  in  the  head-plate.  Since  the  groove  in  the  platform  of  the 
imperfection  measuring  device  is  concentric  with  the  vertical  edge  of  the 
circular  groove  in  the  base-plate,  the  corresponding  edge  of  the  circular 
groove  in  the  head-plate  will  be  concentric  with  its  counterpart  in  the 
base-plate  when  the  pointer  of  the  dial  indicator  is  undisturbed  as  the 
indicator  mechanism  is  moved  smoothly  along  the  platform  groove.  Thus, 
appropriate  alignment  of  the  head-plate  and  the  base-plate  of  the  testing 
fixture  was  accomplished. 
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(11-6).  SPECIMEN  INSTALLATION.  The  procedure  used  to  install  specimens  in 
the  testing  fixture  was  similar  for  specimens  with  unsupported  straight 
edges  and  for  specimens  with  simply-supported  straight  edges.  It  is  con¬ 
venient  to  describe  the  installation  procedure  for  the  specimens  with 
unsupported  straight  edges  first,  and,  subsequently,  describe  the  additional 
installation  features  associated  with  the  specimens  for  which  the  straight 
edges  were  simply-supported.  In  either  case  the  curved  edges  of  all  speci¬ 
mens  were  clamped. 

In  preparation  for  buckling  tests  of  specimens  with  unsupported  straight 
edges,  a  specimen  was  centered  in  the  circular  groove  of  the  base-plate  of 
the  testing  fixture  and  the  pressure  blocks  were  adjusted  by  finger  tightening 
appropriate  screws.  The  platten  of  the  testing  machine  was  then  raised  to 
allow  the  upper  curved  edge  of  the  specimen  to  enter  the  circular  groove  of 
the  head-plate  to  approx ima te 1 y  three-qua r ters  of  the  groove  depth.  Pressure 
blocks  in  the  head-plate  were  then  adjusted  by  finger  tightening  appropriate 
screws  in  the  head-plate.  Generators  at  several  locations  along  the  circum¬ 
ference  of  a  specimen  were  aligned  vertically  with  a  precision  square.  The 
platten  of  the  testing  machine  was  then  raised  until  the  specimen  experienced 
a  compressive  force  of  25  lb.  Pressure  blocks  in  the  base-plate  and  in  the 
head-plate  were  adjusted  to  their  final  positions  by  applying  a  40  in-lb 
torque  to  appropriate  screws  while  the  25  lb  pre-load  was  maintained. 

Preparation  for  buckling  tests  of  specimens  with  a  simply-supported 
straight  edges  was  similar  to  that  for  specimens  with  unsupported  straight 
edges.  A  specimen  was  centered  in  the  circular  groove  of  the  base-plate 
and  its  straight  edges  were  inserted  in  the  slots  of  the  vertical  edge 
supports  (hence  forth  referred  to  as  book-ends)  as  shown  in  Figure  7- 
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Pressure  blocks  in  the  base-plate  and  at  each  book-end  were  adjusted  to 
the  finger  tight  position.  The  platten  of  the  testing  machine  was  then 
raised  to  allow  the  upper  curved  edge  of  the  specimen  to  enter  the  circular 
groove  of  the  head-plate  to  approximately  three-quarters  of  the  groove 
depth.  Pressure  blocks  in  the  head-plate  were  then  adjusted  to  the  finger 


tight  position.  Generators 
at  several  locations  along 
the  circumference  of  the 
specimen  were  aligned 
vertically  with  a  preci¬ 
sion  square.  The  platten 
of  the  testing  machine 
was  then  raised  until  the 
specimen  experienced  a 
compressive  force  of  50  lb. 
Pressure  blocks  in  the 


base-plate  and  in  the  head-  FIGURE  7.  Test  specimen  installed  in  the  base¬ 
plate  and  book-ends  that  supply 

plate  were  adjusted  to  the  s i mp 1 y-suppor ted  edge  conditions. 

AO  in-lb  position,  and  the 


vertical  pressure  blocks  in  the  book-ends  (see  Figure  8)  were  adjusted  to 


the  finger  tight  position  for  specimens  of  Program  A. 

For  Program  B  the  pressure  exerted  on  the  specimen  by  the  vertical 


pressure  blocks  was  adjusted  to  different  amounts. 
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FIGURE  8.  Simply-supported  edge  of  a  test  specimen 
showing  vertical  pressure  blocks. 
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SECTION  I  I  I 


EXPERIMENTAL  DATA 

(lll-l).  INITIAL  GEOMETRICAL  IMPERFECTIONS.  Initial  deviations  of  a 
specimen  from  a  perfect  cylindrical  surface  are  referred  to  in  this 
report  as  initial  geometric  imperfections.  Imperfection  measurements 
were  obtained  at  the  nodes  of  a  rectangular  grid  drawn  on  the  inner  sur¬ 
face  of  a  specimen.  Details  of  the  imperfection  grids  for  specimens 
with  different  aspect  ratios  are  listed  in  Appendices  A  and  B  for  speci¬ 
mens  with  simply-supported  straight  edges  and  unsupported  straight  edges, 
respectively.  This  data  corresponds  to  testing  Program  A.  Details  of 
the  imperfection  grids  used  in  testing  Program  B  are  listed  in  Appendix  C. 

For  testing  Program  A,  the  reference  point  for  imperfection  measure¬ 
ments  lies  on  the  specimen  centerline  two  inches  below  the  upper  clamped 
edge.  Measurements  were  made  at  equal  intervals  on  either  side  of  the 
specimen  centerline,  and  at  equal  intervals  along  the  generators  of  a 
specimen.  Measurements  for  specimens  with  unsupported  edges  were  obtained 
along  these  edges,  while  measurements  for  specimens  with  simply-supported 
edges  were  obtained  as  close  to  the  straight  edges  as  the  imperfection 
measuring  device  would  permit.  These  observations  are  ref'ected  in  the 
imperfection  data  contained  in  Appendices  A  and  E.  Actual  spacing  dimen¬ 
sions  for  the  imperfection  grid  associated  with  each  panel  aspect  ratio 
are  shown  in  these  appendices  also. 

For  testing  Program  B  a  specially  constructed  device,  with  two  circu¬ 
lar  members  that  were  concentric  with  the  groove  in  base-plate,  was  used 
to  establish  the  reference  point  for  imperfection  measurements.  Measure¬ 
ments  for  imperfections  were  made  at  the  locations  described  for  Program  A. 
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Actual  spacing  dimensions  for  imperfection  grids  associated  with  panels 
with  unsupported  straight  edges  and  simply-supported  straight  edges  are 
indicated  in  Appendix  C.  The  imperfection  measurements  for  Program  B  are 
also  presented  in  Appendix  C. 

(111-2).  STRAIN  MEASUREMENTS.  It  was  of  interest  to  determine  if  the 
axial  compressive  end-load  was  applied  uniformly  along  the  curved  edges 
of  a  specimen. 

Program  A.  For  testing  Program  A,  axial  strains  were  measured  along  an 
arc  lying  one  inch  below  the  upper  clamped  edge  on  the  specimen  centerline 
and  at  two  other  equally  spaced  points  on  both  sides  of  the  centerline. 
Electrical  resistance  strain  gages  were  bonded  at  identical  locations  on 
both  sides  of  a  specimen  so  that,  by  means  of  an  appropriate  four-arm 
bridge,  only  axial  strains  were  sensed.  Locations  of  the  strain  gages  on 
the  surface  nearest  the  center  of  curvature  of  a  specimen  can  be  observed 
in  Figure  7-  Appendices  A  and  B  contain  details  for  strain  gage  locations 
for  Program  A. 

Program  B.  For  testing  Program  B,  axial  strains  were  measured  a*-  t .*>•,»  five 
locations  described  in  Program  A  for  specimens  with  unsupported  edges. 

Axial  strains  were  measured  at  only  three  locations  for  soecimens  with 
simply-supported  edges:  on  the  centerline  and  at  two  points  symmetrically 
located  relative  to  the  centerline.  Appendix  C  contains  details  for  strain 
gage  locations  for  Program  B. 

Strain  gage  readings  were  obtained  using  a  Vi shay/El 1 i s-20  digital 
strain  indicator  and  a  Vishay-El 1 i s-21  ten  channel  switching  and  balancing 
unit.  This  equipment  is  shown  in  Figure  9.  Strain  data  associated  with 


viously,  specimens  were  not  perfectly  circular  cylindrical  as  they  emerged 
from  the  mold.  Consequently,  installation  in  the  testing  fixture  caused 
regions  near  the  ends  to  assume  a  circular  cylindrical  shape  of  radius  12 
inches,  while  cross  sections  away  from  the  ends  assumed  noncircular  shapes 
with  the  unsupported  edges  "dishing"  toward  the  center  of  curvature.  Gen¬ 
erally,  the  centerline  and  generators  in  the  regions  on  either  side  of  the 


centerline  were  straighter  than  generators  near  the  edges.  Since  the  axial 
stiffness  of  a  small  strip  of  specimen  parallel  to  a  generator  depends  on 
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its  flexural  stiffness  as  well  as  its  in-plane  stiffness,  the  strains  asso¬ 


ciated  with  strips  located  at  various  sites  on  a  specimen  should  be  expected 
to  be  different.  Indeed,  the  non-uniformity  of  the  strain  distribution 
should  be  expected  to  be  intensified  as  the  loading  process  proceeds.  It  is 
noted,  from  the  strain  data  for  specimens  with  unsupported  straight  edges  in 
Appendix  B,  that  the  strains  near  the  straight  edges  actually  changed  from 
compression  to  tension  during  the  early  stages  of  loading.  This  reversal  of 
strain  near  the  straight  edges  is  associated  with  the  bifurcation  load  for 
the  specimen  under  an  essentially  uniform  end- 1 oad  (as  opposed  to  uniform 
end-displacement).  A  more  detailed  discussion  of  this  behavior  is  presented 
in  a  later  section. 

Now  consider  specimens  with  s impl y-supported  straight  edges.  These 
specimens  experienced  the  same  "dishing"  effects  as  specimens  with  unsupported 
straight  edges;  however,  the  book-ends  forced  the  edges  to  become  straight 
with  the  consequent  configuration  change  on  the  interior  of  a  specimen. 
Accordingly,  since  these  initial  installed  configurations  occur  in  an  essen¬ 
tially  random  manner,  nonuniform  strain  distributions  should  be  expected 
during  the  loading  process. 

Another  possible  reason  for  the  nonuniform  strain  distributions  exhibited 
by  specimens  with  simply-supported  edges  is  the  asymmetry  in  the  axial  dis¬ 
placements  that  can  be  introduced  at  the  book-ends.  The  strain  distribution 
on  the  interior  of  a  specimen  will  deviate  from  uniformity  if  the  axial  dis¬ 
placement  distributions  along  the  two  simply-supported  edges  differ. 

The  further  observation  should  be  made  that  not  only  do  the  initial  imper¬ 
fections  influence  the  buckling  resistance  of  a  specimen,  but  a  certain  unknown 
initial  stress  distribution  is  associated  with  the  installed  initial  configura¬ 
tion  that  also  influences  the  buckling  resistance.  The  influence  of  the 
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initial  stress  distribution  on  the  buckling  resistance  is  rather  capricious. 
The  initial  geometric  imperfections  that  are  incorporated  in  CLAPP  assume 
that  the  imperfect  configuration  is  stress  free. 

(I  I  1-3) •  END-SHORTENING.  The  end-shortening  of  each  specimen  in  testing 
Programs  A  and  B  was  measured  with  a  dial  indicator  that  was  positioned  so 
as  to  determine  the  relative  displacement  between  the  platten  and  the 
cross-head  of  the  Tinius-Olsen  testing  machine.  The  axial  compressive  load 
applied  to  the  specimen  was  read  directly  from  the  testing  machine.  Load 
and  corresponding  end-shortening  data  are  listed  in  Appendices  A  and  B  for 
specimens  of  testing  Program  A  and  in  Appendix  C  for  specimens  of  testing 
Program  B. 

(lll-A).  BUCKLING  BEHAVIORS.  Experimental  buckling  loads  for  the  AO  speci¬ 
mens  associated  with  testing  Program  A  are  listed  in  column  7  of  Table  1 
and  those  associated  with  the  10  specimens  of  testing  Program  B  are  listed 
in  column  5  of  Table  2.  The  recorded  experimental  buckling  load  for  each 
specimen  of  either  testing  program  was  characterized  by  a  distinct  loss  in 
load  carrying  capacity.  Loss  of  load  carrying  capacity  was  detected  readily 
from  the  load-dial  of  the  Tinius-Olsen  hydraulic  testing  machine,  and  was 
always  accompanied  by  a  sudden  shift  in  the  equilibrium  configuration  of 
the  specimen  that  was  clearly  audible  and  visible. 
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SECTION  IV 


EXPERIMENTAL  AND  NUMERICAL  RESULTS 

( I V- 1 )  .  ANALYSIS  OF  RESULTS  OF  TESTING  PROGRAM  A.  The  theoretical  buckling 
loads  shown  in  column  8  of  Table  1  were  calculated  assuming  that  each  panel 
was  subjected  to  a  force  that  was  uniformly  applied  along  its  curved  edge. 

The  bifurcation  branch  of  CLAPP  was  used  to  calculate  these  buckling  loads. 
Consequently,  the  theoretical  buckling  loads  listed  in  Table  1  correspond 
to  bifurcation  under  a  uniformly  applied  axial  compression  using  both  a  12 
x  12  and  a  20  x  12  finite  difference  grid.  For  this  bifurcation  analysis  the 
20  x  12  grid  yielded  a  bifurcation  load  less  than  three  percent  smaller  than 
the  12  x  12  grid.  Nevertheless,  numerical  calculations  for  the  bifurcation 
load  for  each  of  the  remaining  specimens  were  obtained  using  either  a  20  x 
12  or  a  16  x  12  finite  difference  grid.  The  clamped  boundary  condition  along 
a  curved  edge  of  a  specimen  requires  the  transverse  displacements  be  zero 
along  a  line  of  nodal  points  coincident  with  the  curved  boundary  and  along 
a  parallel  line  of  nodal  points  just  inside  the  boundary.  The  larger  num¬ 
ber  of  grid  points  was  always  taken  along  the  generators  of  a  specimen  to 
minimize  this  internal  constraint. 

SPECIMENS  WITH  S IMPLY-SUPPORTEO  STRAIGHT  EDGES.  Table  1  shows  that  the  bifur¬ 
cation  load  predicted  by  CLAPP  was  greater  than  its  corresponding  set  of 
experimental  buckling  loads  except  for  one  specimen.  The  ratio  of  the  modi 
fied  experimental  load  to  the  numerical  bifurcation  load  is  shown  in  the  last 
column  of  Table  1  for  each  specimen.  The  modified  experimental  buckling  load 
was  obtained  by  assuming  that  the  observed  experimental  buckling  load  was 
uniformly  distributed  along  the  curved  edge  of  a  specimen  and  only  the  portion 
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of  the  curved  edge  between  the  book-ends  contributed  to  the  buckling  of  a 


specimen. 

The  ratio  of  the  modified  experimental  buckling  load  to  the  theoretical 
bifurcation  load  fell  in  the  range  0.5  <  p  <  1.0  for  eighteen  of  the  twenty 
specimens  with  simply-supported  straight  edges.  The  ratio  was  0.458  for 
one  specimen  and  1.072  for  another.  We  remark  that  theoretical  bifurcation 
loads  predicted  by  the  energy  method  represent  upper  bounds  to  the  classical 
bifurcation  loads  associated  with  the  test  specimens.  Moreover,  in  the 
presence  of  initial  geometric  imperfections  the  experimental  buckling  loads 
can  be  expected  to  be  less  than  the  classical  bifurcation  load  if  the  speci- 
ment  is  imperfection  sensitive.  If  a  specimen  is  not  sensitive  to  initial 
imperfections  its  experimental  buckling  load  can  be  expected  to  compare 
favorably  with  the  theoretical  bifurcation  load. 

Another  complicating  factor  is  the  nature  of  the  simply-supported  boun¬ 
dary  condition.  Actually  two  different  types  of  the  simply-supported  boun¬ 
dary  conditions  must  be  recognized.  Normally,  a  simply-supported  boundary 
condition  implies  the  bending  moment  and  transverse  displacement  are  zero. 

In  addition  to  these  conditions,  in-plane  displacements  or  membrane  forces 
must  be  specified.  Consequently,  with  reference  to  the  simply-supported 
test  specimens,  circumferential  displacement  can  be  prevented,  can  be 
allowed  to  occur  freely,  or,  there  can  be  some  intermediate  partial  restric¬ 
tion  of  circumferential  displacements.  The  theoretical  bifurcation  loads 
listed  in  Table  1  correspond  to  freely  occurring  circumferential  displace¬ 
ments  along  the  straight  edges. 

The  foregoing  observations  are  offered  in  explanation  of  the  rather 
large  difference  between  the  experimental  buckling  load  and  the  theoretical 
bifurcation  load  for  some  test  specimens  and  the  much  better  agreement 
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between  the  two  loads  for  other  specimens.  The  ratio  1.072  for  one  test 
specimen  is  believed  to  be  a  quirk  arising  out  of  the  way  the  modified 
experimental  buckling  load  is  defined. 

Figure  10  shows  the  experimental  and  numerical  load  versus  end¬ 
shortening  curves  for  a  16  in.  x  16  in.  panel  with  simply-supported 
straight  edges.  The  numerical  load  vs  end-shortening  curve  was  obtained 
using  measured  initial  imperfections.  A  plot  of  the  normalized  buckling 
determinant  is  also  shown  in  the  figure.  We  note  that  the  numerical  curve 
is  essentially  linear  even  though  initial  imperfections  are  present. 
Furthermore,  the  panel  buckled  apparently  by  bifurcation  as  indicated  by 
the  system  buckling  determinant  becoming  negative. 

CLAPP  was  modified  to  provide  the  capability  of  prescribing  uniform 
end  displacements  as  opposed  to  prescribing  uniform  end  load. 

Figure  11  shows  the  experimental  and  theoretical  load  versus  end¬ 
shortening  curves  for  the  16  in.  x  8  in.  specimen  with  the  [0/90],^  fiber 
pattern.  The  theoretical  bifurcation  loads  for  conceptual  models  with 
simply-supported  edges  are  2522  lb  when  c i rcumferen t i a  1  displacements  occur 
freely,  and  4311  lb  when  circumferential  displacements  are  prevented  at  the 
straight  edges.  The  experimental  buckling  loads  for  the  two  test  specimens 
are  2500  lb  and  3240  lb. 

There  is,  perhaps,  better  agreement  between  the  theoretical  bifurcation 
loads  and  the  experimental  buckling  loads  than  is  immediately  obvious.  Con 
sider  the  experimental  curve  labeled  1  in  Figure  11.  The  urusual  shape  of 
this  load  versus  end-shortening  curve  appears  to  be  the  result  of  the  panel 
slipping  circumferentially  in  the  book-ends  at  a  load  near  300  lb.  In  this 
test  it  becomes  clear  that  the  friction  forces  exerted  by  the  vertical 
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pressure  blocks  on  the  specimen  was  not  sufficient  to  prevent  circumferen¬ 
tial  displacements.  Consequently,  at  approximately  300  lb  the  book-ends 
allowed  circumferential  displacements  to  commence.  However,  the  specimen 
edges  apparently  made  contact  with  back  surfaces  of  the  book-ends  and 
c i rcumf eren t i a  1  displacements  were  prevented  from  approximately  1500  lb 
onward.  One  expects  the  experimental  buckling  load  for  this  specimen  to 
fall  between  the  two  theoretical  bifurcation  loads. 

The  behavior  exhibited  by  the  experimental  curve  labeled  2  in  Figure 
II  can  be  explained  in  a  similar  manner.  Initially  the  friction  forces 
exerted  by  the  vertical  pressure  blocks  on  the  specimen  is  sufficient  to 
prevent  circumferential  displacements.  At  a  load  of  approximately  800  lb 
a  gradual  slippage  occurs  and,  finally  at  a  load  of  approximately  2250  lb 
the  friction  forces  were  not  sufficient  to  prevent  or  restrict  these  dis¬ 
placements  any  longer.  A  sudden  slippage  occurred  that  resulted  in  the 
buckling  of  the  specimen  before  its  straight  edges  made  contact  with  the 
back  surfaces  of  the  book-ends.  Accordingly,  one  expects  the  experimental 
buckling  load  for  this  specimen  to  agree  closely  with  the  theoretical 
bifurcation  load  corresponding  to  the  simple-support  condition  when  circum¬ 
ferential  displacements  occur  freely. 

The  experimental  strain  distribution  shown  in  Figure  12  shows  an  en¬ 
couraging  resemblance  to  the  calculated  nodal  load  distribution  shown  in 
Figure  13- 

SPECIMENS  WITH  UNSUPPORTED  STRAIGHT  EDGES.  Table  1  shows  that  the  bifurca¬ 
tion  load  predicted  by  CLAPP  was  less  than  its  corresponding  set  of  experi¬ 
mental  loads  for  every  specimen.  This  unexpected  behavior  can  be  explained 
as  follows.  Because  the  radius  of  a  specimen  was  less  than  12.0  in. 
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distortions  of  the  straight  edges  of  a  specimen  occurred  upon  clamping 
its  curved  edges  in  the  testing  fixture.  The  experiment  buckling  load  is 
believed  to  correspond  to  a  limit  point  load  that  is  greater  than  the 
bifurcation  load  for  a  perfectly  cylindrical  geometry.  In  other  words, 
the  initial  distortion  of  the  straight  edges  of  a  specimen  prevents  an 
experimental  detection  of  a  bifurcation  load  corresponding  to  a  wrinkling 
of  these  edges  for  a  perfect  cylindrical  geometry  under  a  uniformly  dis¬ 
tributed  edge  compression. 

We  remark  that  the  bifurcation  loads  predicted  by  CLAPP  correspond  to 
cylindrical  specimens  without  initial  imperfections  under  forces  that  are 
uniformly  distributed  along  the  curved  edges.  It  is  important  to  note 
that  axial  and  circumferential  displacements  along  the  curved  edges  are 
not  specified.  Consequently,  theory  predicts  a  specimen  will  collapse  at 
the  bifurcation  load. 

The  experimental  model  is  clamped  tightly  in  the  head-plate  and  in  the 
base-plate  so  that  c i rcumferent i a  1  displacements  along  the  curved  edges  are 
prevented.  More  importantly,  axial  deformations  are  restricted  by  the 
relative  movement  of  ihe  cross-head  and  the  platten  of  the  testinq  machine. 
Consequently,  when  the  straight  edges  buckle  the  testing  fixture  prevents 
the  collapse  that  theory  predicts  for  the  specimen  under  uniform  load.  As 
a  result  of  the  clamp  conditions  the  specimen  retains  significant  load 
carrying  capacity.  In  fact,  the  experimental  buckling  loads  correspond  to 
limit  points  in  the  post  buckled  region. 

A  close  inspection  of  the  strain  data  indicates  the  strain  at  the 
straight  edges  changed  sign  at  loads  that  are  in  the  appropriate  neighbor¬ 
hoods  of  the  theoretical  bifurcation  loads.  The  sign  change  on  the  edge 
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strains  corresponds  to  the  observed  severe  distortion  of  the  straight  edges 
of  a  specimen.  The  load  levels  at  which  the  sign  change  occurred  depended 
on  the  initial  imperfection  in  the  specimen. 

CLAPP  was  modified  so  that  the  buckling  behavior  of  specimens  under 
prescribed  uniform  end-displacements  could  be  studied. 

Figure  14  shows  the  experimental  equilibrium  paths  for  the  16  in.  x 
8  in.  specimens  with  the  [ 0/90 J 2 s  f'ber  pattern.  It  also  shows  equilibrium 
paths  predicted  by  CLAPP  for  initial  imperfections  of  the  form  WQ(x,y)  = 

Wl ( l+cosTrx/2a)  for  amplitudes  Wl  =  0.0,  +0.005,  and  -0.010.  The  limit  point 
buckling  loads  discussed  previously  are  clearly  evident.  The  ratios  of  the 
experimental  buckling  loads  to  the  theoretical  limit  load  for  a  perfectly 
cylindrical  specimen  are  0 . 85  and  0.83-  The  ratios  of  the  experimental 
buckling  loads  to  the  theoretical  limit  load  for  an  initial  imperfection 
of  amplitude  0.005  are  0.93  and  0.9'.  Even  better  agreement  is  expected 
for  the  exact  imperfection  distributions. 

Figure  15  shows  the  strain  distribution  at  various  load  levels  for  the 
specimen  discussed  in  the  previous  paragraph.  Evidently,  the  force  applied 
to  the  curved  edges  of  the  specimen  was  uniform  during  the  early  stages  of 
the  loading  process;  that  is,  during  the  initial  increments  of  prescribed 
end-displacements.  Accordingly,  the  bifurcation  load  predicted  by  CLAPP 
for  a  uniformly  applied  end-load  should  agree  qualitatively  with  the  load 
at  which  the  sign  change  occurred  for  the  edge  strains.  The  theoretical 
bifurcation  load  from  fable  1  is  476  lb  and  the  load  at  which  the  axial 
edge  strain  ceased  to  increase  was  near  600  lb. 

Figure  16  shows  the  nodal  forces  predicted  by  CLAPP  for  the  same  speci- 
ment  distribution.  During  the  initial  increments  of  prescribed  end- 
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Analytical  and  experimental  load  vs  end-shortening  curves  for  the  perfect  and  the 
imperfect  panel  with  unsupported  straight  edges. 


experimental  strain 
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displacements  the  strain  distributions  and  the  nodal  force  distributions 
are  nearly  uniform.  These  distributions  become  nonuniform  as  the  loading 
process  continues;  however,  there  is  a  visible  agreement  between  the 
experimental  strain  distribution  and  the  calculated  nodal  load  distribution 
for  all  levels  of  prescribed  end-displacements. 

(IV-2).  ANALYSIS  OF  RESULTS  OF  TESTING  PROGRAM  B.  The  theoretical  buckling 
loads  shown  in  column  6  of  Table  2  were  calculated  assuming  that  each  panel 
was  subjected  to  prescribed  end-displacements  that  were  uniform  along  the 
curved  edges.  Calculations  were  made  for  16  x  11  and  20  x  1 1  finite- 
difference  grids  to  check  the  convergence  of  the  numerical  process. 

SPECIMENS  WITH  SIMPLY-SUPPORTED  EDGES.  The  experimental  equilibrium  paths 
for  five  16  in.  x  12.83  in.  specimens  with  the  [0/+^S/30] ^  fiber  pattern 
and  simply-supported  straight  edges  are  shown  in  Figure  17.  The  load  at 
which  each  specimen  buckled  is  marked  on  the  graph  beside  the  correspond i ng 
equilibrium  path.  The  theoretical  equilibrium  path  is  also  shown  in  the 
same  figure  for  the  perfect  panel  and  for  an  initial  imperfection  amplitude 
of  -0.010  in.  The  two  paths  essentially  coincide  except  that  their  termina¬ 
tion  points  (buckling  loads)  are  different  (806A  and  6608  lb.  for  the  per¬ 
fect  and  imperfect  panels,  respectively).  These  buckling  loads  correspond 
to  a  16  x  11  finite-difference  grid.  A  20  x  1 1  finite-difference  grid  was 
also  used  to  calculate  the  buckling  loads  for  the  same  perfect  and  imperfect 
panels.  The  equilibrium  paths  in  this  case  were  indistinguishable  from 
those  obtained  using  the  16  x  11  finite-difference  grid  except  that  their 
termination  points  (buckling  loads)  are  somewhat  less  (6950  and  5622  lb. 
for  the  perfect  and  imperfect  panels,  respectively)  than  those  associated 
with  the  16  x  11  finite-difference  grid. 
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Figure  18  shows  theoretical  equilibrium  paths  for  the  perfect  panel  and 
the  imperfect  panel  (Wl  =  -0.010  in.)  where,  instead  of  using  end-displacement 
as  a  measure  of  the  system  displacement,  the  norm  of  the  transverse  displace¬ 
ments  is  used.  These  curves  exhibit  more  clearly  the  nonlinear  character  of 
the  equilibrium  path.  These  curves  correspond  to  the  20  x  1 1  finite- 
difference  grid.  Also  note  the  mathematical  form  of  the  geometric  initial 
imperfection  shown  on  the  same  figure. 

The  theoretical  buckling  load  was  detected  numerically  as  a  change  in 
the  system  buckling  determinant  from  positive  to  negative.  This  indicates 
that  the  simply-supported  panel  experienced  bifurcation  buckling.  A  plot 
of  the  normalized  buckling  determinant  versus  total  end  load  for  an  initial 
imperfection  amplitude  of  -0.010  in.  is  shown  in  the  upper  left  hand  corner 
of  Figure  1 8 . 

The  initial  imperfection  amplitude  corresponds  to  a  radial  deviation 
from  a  perfect  cylinder  of  0.020  inches  at  the  geometric  center  of  the  panel, 
which  is  approximately  equal  to  one  half  the  panel  thickness.  This  is  a 
rather  large  initial  imperfection.  The  magnitudes  of  the  measured  initial 
imperfections  were  frequently  larger  than  0.020;  however,  the  distribution 
of  the  imperfections  were  different.  Generally,  the  largest  imperfections 
occurred  near  the  edges  of  a  specimen  with  somewhat  smaller  ones  occurring 
near  the  centerline. 

The  16  x  11  and  20  x  II  finite-difference  grids  gave  8o64  lb  and  6950 
lb  as  the  bifurcation  loads  for  the  perfect  panel.  That  is,  the  20  x  1 1 
grid  predicted  a  bifurcation  load  13-8%  lower  than  the  16  x  11  grid.  Further 
improvement  in  the  theoretical  buckling  load  is  expected  with  a  further 
refinement  in  the  finite-difference  grid.  CLAPP  is  presently  dimensioned  to 
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Figure  17.  Theoretical  and  experimental  equilibrium  paths  for  simply-supported  test 
specimens  for  testing  Program  B. 
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Analytical  load  versus  deflection  curves  for  the  perfect  and  imperfect  panel  with 
simply  supported  straight  edges  and  unconstrained  circumferential  displacements. 


handle  a  maximum  20  x  12  finite-difference  grid,  so  that  further  refine¬ 
ments  were  not  undertaken  in  this  investigation. 

The  16  x  11  and  20  x  1  I  finite-difference  grids  predicted  bifurcation 
loads  of  6608  and  5622  lb  for  the  imperfect  panel  (Wl  -  -0.010  in.).  Thus, 
the  20  x  1 1  grid  predicted  a  bifurcation  load  lb. 3%  lower  than  the  16  x  11 
grid  for  the  imperfect  panel.  Again,  improvement  in  the  theoretical  buckling 
load  for  the  imperfect  panel  can  be  expected  for  a  more  refined  finite- 
difference  grid. 

Figure  18  shows  that  the  experimental  buckling  loads  are  not  grossly 
misrepresented  by  the  buckling  load  predicted  by  the  20  x  1 1  grid  for  the 

imperfect  panel.  Based  on  the  buckling  load  predicted  by  the  20  x  1  I  grid 

for  the  perfect  panel  (6950  lb)  the  ratios  of  the  experimental  to  the  theore- 
tical  buckling  load  for  the  five  specimens  are  0.831,  0.793,  0.7b),  0.716, 
and  0.668.  The  last  ratio  really  should  not  be  counted  because  this  speci¬ 
men  had  been  tested  previously  with  very  loose  straight  edges.  Based  on 

the  buckling  load  predicted  by  20  x  II  grid  for  imperfect  panel  (5622  lb) 

these  ratios  become  1.03,  0.980,  0.918,  0.885,  and  0 . 8 1 8 . 

The  theoretical  curves,  and  thus  the  theoretical  buckling  loads,  upon 
which  the  foregoing  ratios  are  based  correspond  to  simply-supported  straight 
edges  for  which  circumferential  displacements  are  unrestrained.  It  is 
impossible  to  be  certain  of  the  nature  of  the  boundary  condition  along  the 
straight  edges  of  the  test  specimens.  As  an  example  of  the  uncertainty  of 
the  type  of  boundary  condition  that  existed  along  the  straight  edges  of  a 
test  specimen  consider  the  test  specimen  for  which  the  experimental  buckling 
load  was  found  to  be  5775  lb,  the  largest  of  the  bucklinq  loads  of  the  five 
specimens.  For  this  test  the  shoulders  of  the  vertical  pressure  blocks 


were  filed  to  a  rounded  conf igurat ion  so  that  the  blocks  contacted  the  speci¬ 
men  along  a  straight  line.  In  the  other  tests  these  contacting  surfaces 
were  flat  and  tended  to  bend  the  panel  when  the  pressure  was  applied.  The 
former,  it  is  expected,  allowed  a  more  freely  occurring  rotation  at  the 
edge  for  which  the  corresponding  buckling  load  should  be  expected  to  be 
less  than  that  for  the  latter  case.  As  can  be  seen  this  was  not  the  case. 

All  this  suggests  that  the  edge  conditions  at  the  straight  edges  of  the 
test  panel  are  uncertain. 

It  is  felt  that  the  conditions  along  the  straight  edges  of  the  test 
specimens  is  the  principal  source  of  the  difference  exhibited  between  the 
experimental  and  theoretical  buckling  loads.  It  is  also  felt  the  CLAPP  is 
much  less  at  fault  for  the  difference  referred  to. 

It  should  also  be  noted  that  any  pressure  exerted  by  the  vertical  pres¬ 
sure  blocks  on  the  panel  tends  to  retard  the  free  occurrence  of  axial  dis¬ 
placements  along  the  straight  edges.  This  could  give  rise  to  shearing 
stresses  in  the  panel  which  would  lead  to  lower  axial  buckling  loads. 

Finally,  the  clamping  mechan i sm  prevents  circumferential  displacements  from 
occurring  freely  along  the  curved  edges.  This  could  cause  local  distortions 
near  the  curved  edges  that  could  lead  to  lower  buckling  loads.  The  model 
for  which  the  theoretical  buckling  loads  were  computed  assumed  that  axial 
displacements  along  the  straight  edges  and  c i rcumferent ial  displacements  along 
the  curved  edges  are  unimpeded. 

SPECIMENS  WITH  UNSUPPORTED  STRAIGHT  EDGES.  The  experimental  equilibrium 
paths  for  six  16  in.  x  12.83  in.  specimei.  with  the  [ 0 /90 ] 2 s  f'ber  pattern 
and  unsupported  straight  edges  are  shown  in  Figure  19.  The  load  at  which 
each  specimen  buckled  is  marked  on  the  graph  besiae  the  corresponding 


equilibrium  path.  The  theoretical  equilibrium  path  is  also  shown  in  the 
same  figure  for  initial  imperfection  amplitudes  of  -0.0025  and  -0.005-  As 
was  the  case  for  specimens  with  simply-supported  straight  edges,  the  theore¬ 
tical  equilibrium  paths  for  the  16  x  11  finite-difference  qrid  essentially 
coincide  for  initial  imperfection  amplitudes  of  -0.0025,  '0.0050,  and  -0.010 
in.  Moreover,  the  buckling  loads  corresponding  to  these  initial  imperfec¬ 
tion  amplitudes  differed  only  slightly,  as  can  be  seen  from  Figure  20. 

The  equilibrium  path  cor  respond  i  ng  to  a  20  x  1)  f.i  n  i  te-d  i  f  ference  grid 
differs  only  slightly  from  that  cor  respond i ng  to  a  16  x  11  grid.  The 
buckling  loads  for  the  )6  x  1)  and  20  x  II  grids  for  an  initial  imperfec¬ 
tion  amplitude  of  -0.01  in.  are  3*4l6  and  3320  lb,  respectively.  In  this 
case  the  refined  f i n i te-d i f ference  grid  did  not  influence  the  buckling  load 
nearly  as  much  as  it  influenced  the  buckling  loads  for  the  simply-supported 
pane  I  . 

The  experimental  equilibrium  paths  depicted  in  Figure  19  indicate  that 
the  stiffness  of  each  test  specimen  changes  noticeably  in  the  load  range 
*400-1000  lb.  This  change  in  stiffness  coincides  also  with  the  observed 
reversal  of  the  axial  strains  at  the  straight  edges  of  the  test  specimens. 
Consequently,  it  is  believed  that  this  change  in  stiffness  signifies  the 
load- level  at  which  the  straight  edges  of  a  specimen  lost  their  capability 
to  resist  greater  axial  loads. 

The  theoretical  equilibrium  paths  shown  in  Figure  20  for  various  initial 
imperfection  amplitudes  reveal  that  a  bifurcation  occurs  at  approximately 
113^  lb  for  a  perfect  panel  This  agrees  well  with  the  experimental  loads 
corresponding  to  the  changes  in  stiffness  of  the  test  specimens  as  shown 
in  Figure  19. 


Using  the  limit  load  (3320  lb)  corresponding  to  an  initial  imperfection 
amplitude  of  -0.010  and  the  20  x  1!  finite-difference  grid,  the  ratios  of 
the  experimental  limit  load  to  the  theoretical  limit  load  are  0.8t8,  0.753, 
0.7^1,  0.693.  and  0.652.  The  theoretical  limit  load  used  to  calculate 
these  ratios  are  based  on  an  axially  symmetric  distribution  of  initial 
imperfections  (W  =  Wl  ( l+cosiu'/a) )  .  The  actual  distributions  of  initial 
imperfections  indicated  a  pronounced  distortion  of  the  straight  edges  of 
the  general  shape  represented  by  the  mathematical  equation  expressed  above, 
that  dimenished  as  the  specimen  centerline  was  approached.  Actually  the 
center  region  of  the  test  specimens  were  relatively  straight  in  comparison 
to  the  straight  edges.  This  observation  is  revealed  by  the  imperfection 
data  listed  in  Appendix  C. 

The  equilibrium  path  for  an  initial  imperfection  distribution  given  by 

Z 

W  =  Wl  (l+cosnf,/a)  (n/b)  is  shown  in  Figure  20.  This  distribution  more 
o 

closely  represents  the  actual  imperfection  distributions.  This  equilibrium 
path  coincided  with  the  other  paths  shown  in  Figure  20.  The  limit  point  load 
for  this  imperfection  was  3A6O  lb.  which  is  essentially  the  same  as  the  limit 
point  loads  for  the  other  initial  imperfections. 
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Figure  19.  Theoretical  and  experimental  equilibrium  paths  for  specimens 
of  Program  8  with  unsupported  straight  edges. 
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SECTION  V 


LONGITUDINAL  STIFFENERS 

(V-l).  INTRODUCTION.  The  compute'  program  CLAPP  was  modified  to  include 
the  effect  of  longitudinal  stiffeners  on  the  buckling  behavior  of  fiber- 
reinforced  panels  by  adding  the  total  potential  energy  of  the  stiffeners 
to  the  total  potential  energy  of  the  panel.  Accordingly,  the  total  poten¬ 
tial  energy  of  a  stiffener  element  is  developed  in  this  section. 

The  stiffeners  considered  in  this  section  have  thin-wall  open  cross 
sections  and  are  assumed  to  coincide  with  the  finite  difference  grid  lines 

f 

that  are  parallel  to  the  generators  of  the  panel.  It  is  not  necessary 
that  a  stiffener  be  associated  with  every  finite  difference  grid  line. 

A  stiffener  is  assumed  to  be  made  from  an  homogeneous  isotropic  mater¬ 
ial;  however,  the  developments  remain  valid  for  quas i - i sot ropi c,  fiber- 
reinforced  stiffeners.  The  special  properties  of  quas i- i sot rop i c ,  fiber- 
reinforced  stiffeners  are  described  later.  It  is  further  assumed  that 
the  cross  section  of  a  stiffener  does  not  vary  along  its  length,  and  that 
each  stiffener  is  free  of  externally  applied  forces. 

Each  stiffener  is  assumed  to  be  rigidly  attached  to  the  panel  along  a 
finite-difference  grid  line.  A  point  of  attachment  is  defined  as  the 
point  on  the  reference  surface  of  the  panel  that  lies  on  the  normal  through 
the  centroid  of  the  stiffener  cross  section.  Mathematically,  the  displace¬ 
ment  components  associated  with  the  panel  and  the  displacement  components 
associated  with  the  stiffener  are  required  to  be  continuous  along  this 
con  tact  line. 
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(V-2) .  STIFFENER  STRAIN  ENERGY.  An  expression  for  the  strain  energy  asso¬ 
ciated  with  small  displacements  of  straight  beams  with  thin-wall,  open  cross 
sections  is  given  by  Bleich  and  Bleich  [4] 

U  =  2  '  f  El  .  (uV  +  El rr(Vs)*"  +  J +  EA(wcf  +  r(e")  [  dz  (V-l) 

The  quantities  appearing  in  Eq .  (V-l)  are  defined  as  follows.  The  rectangu¬ 
lar  coordinates  t,  and  n  coincide  with  the  principal  centroidal  axes  of 
inertia,  and  z  signifies  a  coordinate  measured  along  the  centroidal  axis  of 

the  beam.  Thus,  l__  and  I  are  principal  centroidal  moments  of  inertia, 

nn 

A  denotes  the  cross  sectional  area,  J  is  the  torsional  constant,  and  r  is 
the  warping  coefficient  for  the  cross  section.  Young's  modulus  and  the 
shearing  modulus  are  denoted  by  E  and  G,  respectively.  Finally,  and  V^ 

are  components  of  displacement  of  the  shear  center  parallel  to  the  E  and  n 
axes,  respectively,  is  the  axial  displacement  of  the  centroid  of  the 
cross  section,  andyfl  is  the  unit  angle  of  twist  that  the  section  experiences. 

The  first  two  terms  in  Eq .  (V-l)  represent  strain  energy  due  to  bending 
about  the  principal  axes  of  inertia,  while  the  next  two  terms  represent 
strain  energy  caused  by  twisting  and  axial  deformation,  respectively.  The 
last  term  represents  strain  energy  associated  with  warping  of  the  cross 
section.  The  strain  energy  associated  with  warping  is  usually  discarded  in 
stiffener  analyses;  it  is  also  discarded  in  the  present  analysis. 

Since  it  is  assumed  that  every  stiffener  is  free  of  externally  applied 
forces,  Eq.  (V-l)  also  represents  the  total  potential  energy  of  a  stiffener. 

Energy  methods  have  been  employed  by  several  investigators  to  examine 
the  effect  of  stiffeners  on  the  buckling  behavior  of  plates  and  panels. 


These  investigators  do  not  agree  universally  on  the  terms  that  need  to  be 
retained  in  the  stiffener  energy  expression.  Donnell  [5]  discards  the 
strain  energies  associated  with  axial  deformations  and  bending  about  an 
axis  perpendi cular  to  the  reference  surface  of  the  panel.  Thus  Donnell 
assumes  that  the  dominant  actions  of  a  stiffener  are  bending  about  an  axis 
parallel  to  a  circumferential  tangent,  and  twisting.  Szilard  [6]  adopts 
the  same  reasoning,  but  argues  that  strain  energy  due  to  twisting  can  be 
discarded  for  closely  spaced  stiffeners.  Palamarchuk  and  Polyakov  [7] 
retain  the  same  bending  and  twisting  energies,  but  include  terms  that 
represent  the  effects  of  externally  applied  forces  and  initial  geometric 
imperfect i ons . 

The  i nves t i ga t ions  cited  consider  only  stiffeners  with  symmetrical 
cross  sections.  Stiffeners  with  unsymmet r i ca 1  cross  sections  are  treated 
in  the  present  developments. 


(V-3).  STIFFENER  ENERGY  IN  MATRI 

X  FORM 
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Connectivity  of  a  stiffener  to  the  panel  is  accomplished  by  expressing 
these  displacement  variables  in  terms  of  the  displacements  (Up,  Vp)  of  the 


*♦7 


point  of  contact  of  a  stiffener 
with  the  panel.  The  situation 


is  depicted  in  Figure  21. 

It  is  shown  in  Ref.  [4]  that, 
for  small  displacements,  the  dis¬ 
placement  components  for  a  generic 


i'i 


point  in  the  rigid  cross  section 


of  a  stiffener  are 


U  =  Us  +  (n  -  n)B 


V  =  V  -  (r.  -  ',)A 

s  s 


FIGURE  21.  General  thin-wall  open 
cross  section. 


At  the  contact  point  D,  f  and  n  =  n  ,  so  that,  with  the  aid  of  Eqs. 

(V-3), 


‘(ns  '  V 
(f;s  "  f’D* 


Because  plane  sections  remain  plane  (r  =  0)  the  axial  displacement  of 
the  contact  point,  wQ,  is 


w  =  w  -  f  U'  -  nn  V* 
D  c  D  c  D  c 


where  u^  and  v^_  are  centroida!  displacements  along  the  K  and  n  axes  as 
shown  in  Figure  (21),  and  (  )'  indicates  differentiation  with  respect  to 


z.  From  Eq.  (V-3) , 


U  =  U  +  n  6  , 
css 


V  =  V  -  n  B  , 
css 


(V-6) 


so  that  Eq.  (V-5)  becomes 
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Eqs .  (V-i»)  and  (V- 5)  lead  to  the  following  transformation 
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Eq.  (V-8)  permits  the  strain  energy  of  a  stiffener  to  be  expressed  in  terms 
of  the  displacements  experienced  by  the  point  of  contact  of  a  stiffener 
wi  th  the  panel . 


(V-4) .  DISPLACEMENT  CONTINUITY.  To  enforce  the  required  displacement  con¬ 
tinuity  along  the  line  of  contact  of  a  stiffener  with  the  panel,  it  is 
necessary  to  project  the  displacemen-  vector  of  the  contact  point  along  the 
circumferential  tangent  and  along  a  normal  to  the  panel.  According  to 
Figure  21 , 
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where  Ut  and  V^  are  the  circumferential  and  normal  components  of  displace¬ 
ment  of  the  contact  point  D. 
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Accordingly,  the  contact 
right-hand  side  of  fcq.  ( V- S) 


point  displacement  vector  appearing  on  the 
is  expressed  as 
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The  assumption  is  made  that  the  panel  thickness 
the  dimensions  of  a  stiffener  normal  to  the  panel, 
tact  point  D  to  lie  in  the  reference  surface  of  the 
Continuity  of  the  displacements  associated  with 
the  stiffener  and  the  displacements  associated  with 
the  panel  requires  that 


is  small  compared  to 
This  permits  the  con- 
panc  I  . 

the  contact  point  on 
the  contact  point  on 


ut  -  V, 


u  =  w, 

n 


wD  =  U, 


8  =  W 


,Y 


8’  =  W 


3"  =  W 


,yxx 


( V- 1 1 ) 


Here  U,  V,  W  are  the  displacement  components  of  a  point  on  the  reference 
surface  of  the  panel. 

Eqs .  (V- 2),  (V- 8),  (V-10),  and  (V-ll)  yield  the  following  matrix  formu¬ 
lation  for  the  strain  energy  of  a  stiffener: 


X, 

v  =  h  j  [d]T[A]T[B]T[S][B][A][d]  dz. 


(V- 12) 
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where 


[A] 


cosa  s i na  0  0  0 

-si  not  cosa  0  0  0 

0  0  10  0 

0  0  0  i  0 

0  0  0  0  I 


(V- 13) 


[B] 


r 


0 

0 


" (ns  '  nD) 


u. 


0 

0 

1 

0 


(V-lM 


[S]  = 


El  0 

nr 

0  EI,U 

0  0 

0  0 


0  0 

0  0 

JG  0 


EA 


(V- 15) 


and 


[d]  =  [V  ,  W  ,  W  ,  W  ,  U  ] 

,  xx  ,xx  ,xxy  ,xy  ,x 

To  interface  with  the  computer  program  CLAPP,  the  elements  of  the 
stiffener  displacement  vector  [d]  are  rearranged  so  that 

>dlT-  >V  \xx-  U,x„-  %>■  ".x,l 

The  total  potential  energy  of  a  stiffener  becomes 


(V- 1 6a) 


(V- 1 6b) 


V  =  % 


J  ^ 


[PAS2]  [d]  dz 


(V- 1 7) 
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The  matrix  [PAS2)  results  from  carrying  out  the  operations  indicated  in 
Eq .  ( V-  1 2 )  taking  into  account  the  alterations  of  the  matrices  [A],  [ B ] , 
and  [S]  because  of  the  rearrangement  of  the  elements  in  the  [d]  matrix. 

The  elements  of  [ PAS2 ]  are  given  in  Figure  22. 

A  final  transformat  ion  is  required  to  express  the  total  potential 
energy  in  terms  of  the  finite-difference  grid  point  displacements. 

(V-5).  FINITE-DIFFERENCE  CONSIDERATIONS.  The  strain  energy  density  of 
a  stiffener  element  is  assumed  to  be  constant  over  the  length  of  the 
element.  Moreover,  the  strain  energy  density  is  assumed  to  be  equal  to 
its  values  at  the  midpoint  of  the  element  length.  Consequently,  the  total 
potential  energy  for  a  stiffener-element  is  taken  as  the  product  of  the 
strain  energy  density  at  the  midlength  of  the  element  and  the  length  of 
the  element.  Symbolically, 

M  =h  [d]T[PAS2] [d]  x  i  ( V- 1 8) 

e 

where  [d]  is  the  displacement  vector  associated  with  the  midpoint  of  the 
s  t i f fene  r-e 1 emen  t . 

Nine  different  stiffener-elements  are  required.  These  stiffener- 
elements  correspond  to  the  nine  types  of  area-elements  used  in  CLAPP;  one 
interior  stiffener-element,  four  boundary  stiffener-elements,  and  four 
corner  stiffener-elements.  The  nine  stiffener-elements  are  shown  in 
Figure  23,  which  also  shows  the  stiffener-element  orientation  and  midpoint 
along  with  the  local  numbering  system  for  the  surrounding  finite-difference 
grid  points. 

For  each  stiffener-element  the  midpoint  displacement  vector  (Eq.  ( V— 1 6b) ) 
is  transformed  to  the  local  grid  point  displacement  vector 
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( v- 19) 


[q]T  =  [wr  ur  vr  w2,  u2,  v2>  ....  w9,  ug,  vg] 

through  the  relations 

d.  =  c. .  q.  (V-20) 

•  U  J 

where  the  matrix  [C]  is  composed  of  coefficients  that  depend  on  the  vari¬ 
able  spacing  of  the  finite-difference  grid.  Accordingly,  the  total  poten¬ 
tial  energy  for  a  stiffener-element  is 

V  =  >a  [q]T  [C]T  [ PAS2 ]  [C]  [q]  *  t  (V-21) 

The  transformation  matrix  [C]  depends  on  the  location  of  the  stiffener- 
element  in  the  global  system.  To  illustrate  the  procedure  used  to  estab¬ 
lish  these  matrices,  the  transformat  ion  matrix  for  the  interior  stiffener- 
element  shown  in  Figure  2A  is  derived. 

In  the  longitudinal  direction,  the  midpoint  (subsequently  referred  to 
as  the  stiffener-element  centroid)  of  the  stiffener-element  contains  the 
centroid  of  the  corresponding  area-element.  This  follows  from  the  defini¬ 
tions  of  an  area-element  and  of  a  stiffener-element.  The  corners  of  an 
area-element  are  the  centroids  of  the  areas  contained  between  adjacent 
grid  lines.  The  length  of  the  corresponding  stiffener-element  is  equal  to 
the  dimension  of  the  area-element  parallel  to  the  generator  of  the  panel. 

A  derivative  of  a  centroidal  displacement  with  respect  to  the  lonqitu- 
dinal  direction  x,  denoted  by  (  )',  can  be  expressed  as  a  linear  combina¬ 

tion  of  displacements  at  the  three  finite-difference  grid  points  (i-1,  i, 
i+l)  lying  on  the  longitudinal  grid  line. 

A  general  one-dimensional  centroidal  function,  f,  and  its  derivatives 
are  expressed  as  linear  combinations  of  the  function  values  at  the  three 
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node  points  on  the  j th  grid  line  as 


(k-h) (3h+k)  ,  (h+3k) (3h+k)  f  _  (k-h) (3k+h) 

1 6k (h+k)  i  +  ]  ] 6hk  i  1 6h (h+k)  i  —  1 ’ 


(V- 22) 
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and 
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f 11  =  _  L  -r-  f 

i  k  (h+kl  i  +  1 


h  (h+k)  ^  i -  1 , 


(V-2A) 


where  h  and  k  are  longitudinal  spacings  between  the  grid  lines  i-1  and 
i,  and  i  and  i+1,  respectively. 

A  derivative  of  a  centroidal  displacement  with  respect  to  the  circum- 
ferential  coordinate  y,  denoted  by  (  ),  can  also  be  expressed  as  a  linear 

combination  of  di splacements  at  suitable  finite-difference  grid  points  by 
means  of  a  Taylor  series.  Accordingly,  a  Taylor  series  expansion  about 
the  point  (i,  j)  yields  the  two  linear  equations 
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(V-25) 


and 


f.  .  =  f.  -  l  f.  +  h  f.,  ( V—  26 ) 
J-'  J  J  J 

where  ?.  and  m  are  spacings  between  the  j-1  and  j,  and  j  and  j  +  1  circumfer¬ 
ential  grid  lines.  Eqs.  (V-25)  and  ( V- 26 )  lead  to  the  first-order  central 
difference  formulas 
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and 


'fj  ■  fj+i  ■  «§  fr  fj-i  (v-28 

The  centroidal  displacement  vector  [d]  includes  mixed  derivatives  also. 
The  coefficients  in  the  transformation  matrix  [C]  that  correspond  to  deri¬ 
vatives  with  respect  to  x  are  obtained  from  Eqs.  (V-23)  and  (V~2k)  for  an 
interior  stiffener-element.  The  coefficients  in  the  transformation  matrix 
[C]  corresponding  to  the  mixed  derivatives  are  determined  by  appropriate 
combinations  of  Eqs.  (V-23),  (V-2A) ,  and  (V- 26). 

a. 

A  general  two-dimensional  centroidal  function,  g,  has  the  mixed  deri¬ 
vative 


g 


]_  • 

2k  9i  +  l 


+ 


!_  I_ 

2h  ”  2k 


•)  9; 


(V- 29) 


Using  Eg.  (V-26)  yields 
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(V-32) 


Identification  of  the  nine  surrounding  grid  points  with  the  local  num¬ 
ber  system  determines  the  coefficients  in  the  transformation 


g  =  c.  g. 


i  -  1,  2, 


9. 


(V- 33) 
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The  coefficients  for  the  mixed  derivative  g"  are  obtained  in  a  similar 
manner . 

The  total  potential  energy  for  any  stiffener-element  is  expressed 
Vg  =  >*  tq]T  [CAS2  ]  [q]  x  i  (V- 

where 
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t 


[CAS2]  =  [C]T  [PAS2]  [C]  ( V- 35) 

Eq,  (V-3M  is  the  formulation  of  the  total  potential  energy  for  any 
stiffener-element  that  interfaces  consistently  with  the  area-elements  used 
in  CLAPP. 

(V- 6).  INCORPORATION  INTO  CLAPP.  The  incorporation  of  lonciitudinal  stiffe¬ 
ners  in  the  computer  program  CLAPP  is  accomplished  with  three  subroutines. 

The  first  of  these  three  stiffener  subroutines  is  named  SCOEF.  SCOER 
calculates  the  elements  of  the  matrix  [ P AS2 3 .  Since  [PAS2J  depends  only  on 
the  material  and  geometric  properties  of  the  stiffener  cross  section  it  is 
calculated  once  for  each  element  and  stored. 

SCOEF  permits  the  user  to  select  any  of  seven  stiffener  cross  sections. 
It  also  provides  for  any  other  cross  section  through  a  user's  choice  ontion. 
For  the  various  stiffener  cross  sections  included  in  SCOEF  see  the  user's 
manual  in  APPENDIX  D. 

The  second  stiffener  subroutine  is  named  STFDIF.  It  calculate^  the 
elements  of  the  trans forma t ion  matrix  [C]  of  Eg  ( V- 20 ) .  The  subroutine 
FDIFF  of  CLAPP  performs  similar  calculations  for  each  area-element.  Con¬ 
sequently,  when  FDIFF  is  called  to  calculate  [C]  for  an  area-element 
STFDIF  is  subsequently  called  from  FDIFF  to  calculate  [C]  for  the  corres- 


59 


ponding  stiffener  element. 

The  third  stiffener  subroutine  is  named  ATBAS.  It  performs  the  matrix 
mu  1 1 i p 1 i ca t i ons  indicated  in  Eq  ( V- 35 ) -  This  subroutine  is  called  from 
the  subroutine  STFD1F.  ATBAS  has  been  constructed  to  take  advantage  of 
the  numerous  zeros  that  occur  in  the  matrix  [C] . 

The  main  program  and  the  subroutine  FDIFF  of  CLAPP  are  modified  so 
that  longitudinal  stiffeners  can  be  included  when  appropriate.  Additions 
to  the  main  program  read  in  stiffener  indices  indicating  (a)  that  stiffe¬ 
ners  are  to  be  included,  (b)  the  choice  stiffener  cross  section,  (c)  whether 
stiffeners  are  located  on  the  inside  or  outside  of  the  panel,  and  (d)  if 
stiffeners  occur  on  every  f i n i te-d i f ference  grid  line.  The  main  program 
calls  SCOEFF  to  construct  the  matrix  [PAS2] .  The  additions  to  the  subrou¬ 
tine  FDIFF  determine  when  the  stiffener  energy  is  required  in  an  analysis 
that  does  not  include  stiffeners  at  every  grid  line.  When  it  is  appropriate, 
STFDIF  is  called  to  calculate  the  stiffener-element  energy  which  is  imme¬ 
diately  added  to  the  correspoinding  area-element  energy.  The  procedure 
used  by  CLAPP  to  assemble  the  system  stiffness  matrix  for  area-elements 
proceeds  in  precisely  the  same  manner  when  stiffeners  are  present. 

(V-7) .  QUASI- ISOTROPIC  FIBER-REINFORCED  STIFFENERS.  A  strain  enerqy 
expression  for  a  quas i - i sotrop ic ,  f iber-reinforced  stiffener  that  is  ana- 
logus  to  the  strain  energy  expression  for  an  isotropic  stiffener  is  estab¬ 
lished  in  this  section.  For  the  purpose  of  the  present  development  a 
quas i - i sot rop i c  material  is  defined  as  one  possessing  midplane  symmetry 
and  alternate  zero  and  90  degree  fiber  directions.  Thus,  for  a  stiffener, 
the  zero  fiber  direction  is  assumed  to  coincide  with  the  longitudinal  axis 
of  the  stiffener,  and  the  90  degree  fiber  direction  is  perpendicular  to 
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this  axis. 


The  stress-strain  relations  [2]  for  a  general  fiber-reinforced  laminate 


are 
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Formulas  connecting  the  elements  of  the  laminate  stiffness  matrices  [A], 
[B],  and  [D]  with  the  material  properties  and  the  geometric  locations  of  the 
individual  layers  are 


Au  -j.V  (hk- vi>- 


B. .  «  ~  Z  Q. -k  (h*  -  b2  ), 
ij  2  k=)  ij  k  k-i 


(v- 38) 


(V-39) 


D.  .  =  l  Q.  ,k  (h,3  -  h,3  (i,j  =  1,2.6) 


U 


k=l 


'ij  '  k  k-1 


(V-kO) 


Pertinent  geometric  quantities  appearing  in  these  relations  are  defined  in 
Figure  25- 

The  laminate  is  assumed  to  possess  midplane  symmetry  so  that  the  B.^. 
are  identically  zero  in  Eqs .  (V-36)  and  CV- 37)-  It  is  also  assumed  that 


6! 


the  stiffener  experiences  a  membrane  state  of  stress  like  the  one  shown 

in  Figure  26.  Consequently,  the  dominant  bendi  >g  action  results  from  the 

moment  caused  by  the  membrane  force  N  . 

x 

If  the  cross  section  of  a  stiffener  is  assumed  to  be  riqid,  then  e  =0 

y 

so  that 

Nx  =  Ail  ex  +  A1G  (2exy)  (V-M) 

and 

N  =  A16  e  +  A66  (2e  )  (V-4 2) 

xy  x  DD  xy 

Now  Aj6  "  0  because  of  the  laminate  fiber  directions  are  either  parallel 
or  perpendicular  to  the  longitudinal  axis  of  the  stiffener.  Consequently, 
the  stress-strain  relations  appropriate  to  an  analysis  of  f i ber- re i nforced 
stiffeners  under  the  foregoing  restrictions  are 

Nx  =  Auex  (V-4  3) 

and 

N  =  Ac6  (2e  )  .  (V-M*) 

xy  00  xy 

These  relations  are  analogous  to  the  isotropic  s t ress-s t ra i n  relations 
used  by  Bleich  and  Bleich  [A]  to  arrive  at  their  expression  for  the  strain 
energy  of  stiffeners  made  from  isotropic  materials.  All  that  is  required 
to  make  the  isotropic  strain  energy  expression  valid  for  quas i - i sot rop i c 
stiffeners  is  to  appropriately  identify  E  and  G  of  the  isotropic  case  with 
Ajj  and  A^^  for  the  quas i - i sot rop i c  case. 

To  do  this,  consider  that 
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(V-4 5) 


Ah  =  s  Qnk  (hk  -  hk.,)  =  t  2  Qnk 


k=l 


k=l 


and 


A66  =  Z  Q66k  (hk  -  hk_,)  =  t  I  Q66k 


(V-46) 


k=l 


k=  I 


if  the  thickness  of  each  layer  is  the  same. 

The  subscript  k  signifies  that  the  k~th  layer  is  under  consideration. 

The  Q.j  are  symmetrical  and  are  the  materia!  coefficients  referred  to  a 
generic  set  of  axis.  They  are  expressed  in  terms  of  the  material  coeffi¬ 
cients,  Qjj>  associated  with  the  material  axes  of  a  layer  through  the  trans¬ 
format  ion: 


Ql  i  =  Qucos4  0  +  2  (Q j 2  +  2  Q66)  sin2  0  cos2  0  +  Q?2  sin4  6, 

Q12  =  (Ql  1  +  Q??  ~  A  Qg6 )  sin2  0  cos2  0  +  O42  (sin4  0  +  cos4  6), 

Q22  =  Ql  1  sin4  0  +  2  (Qj  o  +  2  Q66)  sin2  0  cos2  0  +  Q2?  cos4  0,  (V- A  7 ) 

Ql 6  =  (Q11"  Ql 2  ■  2  Qs6)  sin  0  cos3  6  +  (Q12  -  O.22  +  2  Q66)  sin3  6  cos  0, 

Q26  =  (Qu  ■  Ql?  ■  2  Q6 6 )  sin3  0  cos  0  +  (Q12  -  Q22  +  2  Q6g)  sin  6  COS3, 

Q66  =  (Ql 1  +  Q2?  ■  2  Ql?  “  2  Qoe)  sin2  0  cos2  0  +  Q66(sin4  0  +  cos4  6). 


The 

through 


quant i t ies  Q.  . 

'  J 


the  relations: 


are  related  to  the  engineering  material  constants 
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Qn  = 


' 'v12v?  1 


Ql  2  =  Q21  = 


v1?  Eu  ^2  1^2  2 
1  -  v  i  ->  \>2 1  1  -  v  j  2  v  2  ]  * 


(V-I48) 


l-Vj2V2l 


$66  -  G12- 


E11,  ^22  are  Young's  modulii  of  elasticity  parallel  and  perpendicular  to 
the  fiber  direction,  respectively;  Gj2  is  the  shearing  modulus  of  elasti¬ 
city  associated  with  the  directions  parallel  and  perpendicular  to  the  fiber 
direction;  and  v]2,  v2  i  are  Poisson  ratios.  vi2  is  associated  with  a 
strain  in  the  1-direction  due  to  a  stress  in  the  2-direction  and  v2]  is 
associated  with  a  strain  in  the  2-direction  due  to  a  stress  in  the  1- 
direction.  It  is  convenient  to  let  the  1-direction  coincide  with  the 
fiber  di recti  on. 

From  Eqs.  (V- 47 )  and  ( V- 48 ) 


k  tu 

Qi  1  =  t -  for  fibers  parallel  to  the  stiffener  axis 

1  1-V12V21  V 


—  k 

Qn  =* 


(V-i»  9) 


k 

0.2 ?  =  “ —  for  fibers  perpendicular  to  the  stiffener 

axis 


Qf,6  ~  066  "  G12  » 


(V-50) 


for  either  fiber  direction. 
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If  N  is  the  number  of  layers  in  the  laminate,  then  from  Eqs.  (\J-Uk) 


and  (V-49) 


An  - 


Nt 

| " v 12^2  1 


Ell  +  E22 


+  X 


where 


0 


X  -< 


1 


i f  N  is  even 

if  N  is  odd. 


(V-51) 


(V-52) 


Note  that  when  N  is  odd,  the  middle  or  odd  layer  is  assumed  to  be  parallel 
to  the  stiffener  axis.  If  the  middle  or  odd  layer  is  perpendicular  to  the 
stiffener  axis  En  and  E22  interchange  positions  in  Eq  (V- 51). 

Fin.'  i  ly, 


A66  -  Nt  G12 


(V-53) 


for  N  even  or  odd.  Notice  that  Nt  is  the  total  thickness  of  the  laminate 
from  which  the  stiffener  is  constructed. 

Now  if  E  and  G  in  the  isotropic  strain  energy  expression  for  thin-wall 
open  sections  is  replaced  with 


All 

JTT 


Ell 


+  E 
2 


22 


+  X  < 


(V-54) 


and 


G=> 


=  G 


I  2  • 


(Y-55) 


one  arrives  at  a  strain  energy  expression  for  the  quas i - i sot rop i c  material 
considered  in  this  section. 
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APPENDIX  A 

STRAIN,  END-SHORTENING,  AND  IMPERFECTION  MEASUREMENTS  FOR 
TEST  PANELS  WITH  S  IMPL Y-SUPPORTED  STRAIGHT  EDGES 
CLAMPED  CURVED  EDGES 

TESTING  PROGRAM  A 
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etter  in  a  specimen  designator  (BCP-98 1 0-3-3" 1 )  indicates 
e  pattern  [0/+45/90]s  and  the  letter  B  denotes  the  patter 

PANELS  HAD  THE  SAME  DESIGNATION  LABEL. 
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IMPERFECTION  MEASUREMENTS 
(x  10"3  in) 


MENTAL  DATA  FOR  16  x  16  PANELS  WITH  S I MPLY- SUPPORTED  STRAIGHT  EDGES 
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EXPERIMENTAL  DATA  FOR  16  x  8  PANELS  WITH  SIMPLY-SUPPORTED  STRAIGHT  EDGES 
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scontinuity  in  end-shortening  vs. 
could  not  be  readily  distinguished 
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co  o  cr\ 
<a  -4  co 
vo 


la  I  -4  r^-  la  cn  n  in  -  4  -j  \D  o  la  >— 

O  4  4  (A  lD  ^  —  O  ^  -4  <A  CO  1^-  LA 

CMrALALA\£>r^CO<A<AOO  —  N  (A 


0  —  0 
o 

LA  lA  VO 


fA  N  N 
CM  |  CO  I  LA 


vO  N  CO  4  O  vO  O  O 

VO  N  N  CO  <n  (T>  O  - 


o  c  o  o  o 

LA  LA  O  LA  O 

cm  la  r**.  O 


lining 


o  o  oooooooooooo 

o  o  oooooooooooo 

CM  I  CA  I  <T  LA  vO  N  CO  <T\  O  —  CM  co  -4  LA 

• —  — '  —  — 1  •“  CM  CM  CM  CM  CM  CM 


CO  LA  CA  O 

-a-  ca  r^- 

LA  LA  \D 


-a-  —  r*^  \o 

O  M3  CA  /-4 

ca  ca  ca  o 


CAr"*0-4LA-4cMv£>0 
CM  LA  CO  O  CA  O  -4  O 

co  co  oo  cacao 


LA  -4  CM  —  (A  -4  v£>  f".  -4  —  sO  00  —  O  co  <A 

'cm  rA  4  4  CO  vfl  LALACACAvO  0-4  I^CA 

LAvO^OOCOCACAOOO  — 


O  CO  CO  vO  CO  <A  —  CM  O  CA 

■ —  CM  fA  4  vO  CO  03  CO  A- 

r*^  GO  CA  O 


CM  4  vO  LA  4  —  CM  00 

vO  CA  CM  O  vO  0-4 

<—  —  CM  CM  CA  CA  -4  -4 


o  ^  cm  vO  • —  o  r*-  o  o  la  co  o  f^*  ca  o  cm  — 

CM  rAlAr^.COO  O  CA  O  -4  OO  —  LA  CO  CM  LA  CA  CM 

o  —  —  CA  <A  (A  4  *4"  4  LA  LA  LA  v£> 


90 


BCP  992 


IMPERFECTION  MEASUREMENTS 
(x  10"3  in) 


1 


BCP  98) 7-B-7-  1  BS 


BCP  981 7-B-7-2  BS 


ciaiiQii 


APPENDIX  B 

STRAIN,  END-SHORTENING,  AND  IMPERFECTION  MEASUREMENTS  FOR 
TEST  PANELS  WITH  UNSUPPORTED  STRAIGHT  EDGES 
AND  CLAMPED  CURVED  EDGES 


TESTING  PROGRAM  B 
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8 CP  98 !0-8-^-2  C F 


BCP  98I0-B-4-]  CF  (cont'd)  BCP  9810-B-4-2  CF  (cont'd) 
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-3  -A  :+ 1 2  +20 


BCP  98 1 0-B-^t-  1  CF 


BCP  9810-B-4-2  CF 


ED  STRAIG 


3050 


IMPERFECTION  MEASUREMENTS 
(x  10"3  in) 


BCP  9921-A-6-1  EF  BCP  992 1  - A-6-2  EF 


BCP  98l 7-B-5-  I  EF  BCP  98 1 7-B-5-2  EF 


This  imperfection  data  is  questionable.  Imperfection  measuring  device 
was  not  mounted  accurately  enough. 


B-3  EXPERIMENTAL  DATA  FOR  16  x  !6  PANELS  WITH  UNSUPPORTED  STRAIGHT  EDGES 
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MENTAL  DATA  FOR  16  x  12  PANELS  WITH  UNSUPPORTED  STRAIGHT  EDGES 
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These  data  were  taken  before  the  large  aluminum  lock-nuts  were  installed  on  the  Tinius-Olsen  testing  machine. 
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These  data  were  taken  before  the  larqe  aluminum  lock-nuts  were  installed  on  the  Tinius-Olsen  testinq  machine. 


*  This  imperfection  data  is  questionable.  Imperfection  measurinq  device 
was  not  mounted  accurately  enough. 


MENTAL  DATA  FOR  16  x  8  PANELS  WITH  UNSUPPORTED  STRAIGHT  EDGES 
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APPENDIX  C 

STRAIN,  END-SHORTENING,  AND  IMPERFECTION  MEASUREMENTS 
FOR  TEST  PANELS  FOR  TESTING  PROGRAM  B 


FIGURE  Locations  for  strain  and  imperfection  measurements  for  16  x  I? 

panels  with  [0/90]2s  fiber  pattern  and  unsupport°d  straight 
edges.  For  16  x  12  panels  with  [0/+kS/30]s  fiber  pattern  and 
simply-supported  straight  edges  only  the  axial  strains  at  the 
three  interior  locations  were  recorded,  and  the  imperfect  data 
were  recorded  for  the  seven  interior  imperfection  grid  lines. 
(Circumferential  dimensions  shown  correspond  to  the  panels 


Imperfection  measurements  (xlO'^in)  for  16x12  panels  with  fiber  pattern 
[0/+45/90]s  and  simply-supported  straight  edges  (curved  edges  clamped). 

Collapse  load  in  parentheses. 

rounded  edges 

DS-A9-1  (4600  1b)  freely  DS-A9-2  (4975  1b)  30  in- lb  on  side  bars 


r 

OO 

8 

BB 

3 

m 

n 

;  7  7 

3 

B9B9 

3 

1 

CM 

2 

-1 

-2 

2 

BB 

-2 

IBI 

0 

-1 

-2 

2 


0 


5 


3  -3 


-2 


-2 


0 


DS-A10-1  (5775  lb)  10  in-lb 


2 


0 


-2  3 


-6  -2  0 


2  0 


DA- A 1 1-1  (5160  lb)  nearly  freely 


-6 

r 

-7 

-3 

0 

0 

-6 

-2 

-6 

0 

5 

-1 

-8 

-8 

HDD 

10 

4 

-6 

-7 

Ml 

Lz 

10 

3 

-8 

3H 

10  *- 10 


-6  i  -4  -2 


3  3 


-2  -2 


DS-A10-2  (5510  lb)  30  in-lb  on  side  bars 


19 


33  22 


BBBB 

IBI  IB! 

30 

lil 

BBS 

0 

0 

3 

3 

M  M 

30 

5  8  12  9 


2 


12  10 


15 


2  18  22 


-3  lb  17 


0  21 


13 


6  8 


Imperfection  measurements  (xl0~3jn)  for  16-12  panels  with  fiber  pattern 
[0/90]s  and  unsupported  straight  edges  (curved  edges  clamped).  Collapse 
load  in  parentheses. 


DS-BI1-1  (2460  lb) 


8  8 


9 


8  9  11110  111 


8  10  12 


6  8  8 


2 


-7  -2 


6  10 


6 


-3  -3  -3  -3  -3  -3 


DS-B11-2  (2715  lb) 


1 

6 

6 

7 

7 

7 

8 

1! 

H 

3 

6  j  8 

9 

6 

13 

-2 

2 

bh 

9 

8 

_ 

10 

12 

-8  -2 


-12 


-14  -9 


-10 


6 


5  3 


8  |  11 


DS-B10-1  (2300  lb) 


OS- B 1 0-2  (2495  tb) 


8  3  0 


1 1 


8  13  12  13  16 


6  10  18  14 


5  11  8 


9  13  21 


10  I  12  I  11  I  13  17  38  33 


5  6  6  8  13  21 


8  I  I  16 


7  I  10 


DS-B9-1  (2165  lb) 


-6  |  -2  |  3 


DS-B9-2  (2410  lb) 


5 

5 

8 

7 

8 

10 

6 

5 

3 

2 

5 

6 

3 

3 

13  14 


9  12 


6  6 


6  7 


3  6 


10 

8 

6 

2 

-21-2  2  4  81016 


2  4  7  I  10  I  10 


7  I  12  11 


5 


DS-B1I-2  [0/90]„  Straight  Edges  Unsupported  DS-B11-1  1 0/90 ]  Straight  Edges  Unsuppo 


DS- B 11-2  [0/90]_  Straight  Edges  Unsupported  DS-B11-1  [0/90].  Straight  Edges  Unsupported 


DS-810-1  [0/90]  Straight  Edges  Unsupported  DS-B10-2  [0/90]  Straight  Edges  Unsupported 


DS-BIO-1  [0.90],,  Straight  Edges  Unsupported  DS-B10-2  [0/90].  Straight  Edges  Unsupported 


DS-B9-1  [0/90] _  Straight  Edges  Unsupported  DS-B9'2  [0/90 ]  Straight  Edges  Unsupported 


DS-B9-)  [ 0/90 ]  Straight  Edges  Unsupported  DS-B9~2  [0/90]„  Straight  Edges  Unsupported 


Straight  Edges  Straight  Edges 

DS-A9-1  [0/+45/90]  Simply- Supported  DS-A9-2  [0/+45/90]  Simply-Supported 


Cont i nued 

Straight  Edges  Straight  Edges 

DS-A9-I  [0/+45/90]  Simply-Supported  DS-A9-2  [0/+A5/90]s  Simply-Supported 
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Continued 

Straight  Edges  Straight  Edges 

DS-AJO-I  [0/+i»5/90]s  Simply-Supported  DS-A10-?  [0/+^5/90]^  Simply-Supported 


3 


Con t i nued 

Straight  Edges  Straight  Edges 

DS-A11-2  [0/+45/90]  Simply-Supported  DS- A 11-2  [0/+45/90]  Simply-Supported 
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APPENDIX  D 


USER'S  MANUAL  FOR  COMPUTER  PROGRAM  CLAPP 


f 
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CLAPP  USERS  MANUAL 


1.  TITLE  CARD  (18A4) 
Columns  1 - 72 


Problem  description 


2.  CONTROL  CARD  (815) 
Columns  1-5 


6-10 

11-15 


I  RUN  -  Takes  integer  values  1,  2,  3, 

according  as  the  completion  of  an 
analysis  requires  1,  2,  3,  .... 
submissions  of  the  program. 

I  RUN  =  1  for  the  first  submission, 

J  RUN  =  2  for  the  second  submission, 
and  so  forth.  Valid  for  nonlinear 
analysis  only  (IBIF  =  0) .  IRUN 
^ifor  any  analysis  that  is  a  con¬ 
tinuation  of  a  previous  application 
of  CLAPP. 

perfect  panel  or  plate 
imperfect  panel  or  plate 

nonl inear  analysis 
linear  bifurcation  analysis 


I  EX  = 


IBIF  = 


16-20 


NFLAT 


{0  flat  plate  analysis 
1  curved  panel  analysis 


21-25 


NAUTO 


0  x  and  y  coordinates  must  be 
punched  for  input.  Variable 
gr i d  capab i 1 i t ies . 

1  Automatic  generation  of  x  and 
y  coordinates  with  uniform,  but 
possibly  different,  spacings 
for  the  x  and  y  directions. 


0  initial  imperfections  are  repre 
sented  by  a  mathematical 
express  ion . 


I 


26-30 


LGRNG 
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1  Initial  imperfections  are  repre- 
!  sented  by  a  two-dimensional 
j  Lagrange  interpolating  function. 


31-35 


0  When  the  load-deflection  curve 
is  generated  in  seqments  the 
tape  NTAPE  contains  the  infor- 
i  mat  ion  required  by  the  program 
;  to  continue  the  calculations 
j  beginning  with  the  last  load- 
i  level  for  which  a  converged 
|  solution  was  obtained.  To  con- 
!  tinue  calculations  set  INDEX  = 
i  2,3,-.-  and  NUSTRT  -  0. 


NUSTRT  = 


1  It  may  happen  that  it  is  desir¬ 
able  to  adjust  the  loadinq 
sequence  to  begin  calculations 
for  a  second  segment  of  the  load- 
deflection  curve.  To  do  this  set 

INDEX  =  2,  3 . NUSTRT  =  1  and 

reset  the  input  data  for  PZSTRT, 
PXSTRT ,  PYSTRT ,  AND  ZlNCR,  XINCR, 
YINCR. 


36-40 


NCURVE  =  0 


3.  CONVERGENCE  CRITERION  CARD  (F10.0,  215> 

Columns  1-10  EPSI  -  Convergence  criterion  that  determines 

acceptable  displacements  at  a  given 
load-level  (order  of  10  l>) . 

11-15  UMAX  -  Program  terminates  if  convergence  at 

any  load-level  has  not  occurred 
within  UMAX  iterations. 

16-20  LEMAX  -  Program  terminates  after  processing 

LEMAX  load-levels  as  collapse  or 
bifurcation  has  not  occurred. 

Resubmit  program  with  I  RUN  =  2. 


4.  GRID  PARAMETERS  CARD  (815) 


Columns  1-5 


h 


NSCHM  =3 


Finite-difference  grid 
numbered  consecutively 
cross-sect ional  circle 


0  Fi ni te-d i f ference  grid 
numbered  consecutively 
generator 


points 
alonq  a 


points 
alonq  a 
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6-10 


ICTOUT  = 


0  Panel  has  no  cutout 


1  Panel  has  a  cutout 


11-15 

NCOLS 

-  Number 

of 

columns 

of  grid  points 

16-20 

NROWS 

-  Number 

of 

rows  of 

grid  points. 

21-25 

IROWI 

-  Number 

of 

the  row 

in  which  the 

side  of  the  cutout 

nearest  the 

first 

row 

of  grid 

points  lies. 

26-30 

IR0W2 

-  Number 

of 

the  row 

in  which  the 

side  of  the  cutout  furtherest 
from  the  first  row  of  grid  points 
lies. 


3 l - 35  IC0L1  -  Number  of  the  column  in  which  the 

side  of  the  cutout  nearest  the 
first  column  of  grid  points  lies. 

36*^0  IC0L2  -  Number  of  the  column  in  which  the 

side  of  the  cutout  furtherest 
from  the  first  column  of  grid 
points  lies. 


nrows 


1 


I  ROW I  IROW2 


|  ICOL1 

NCOLS 
I  ICOL2 


- f 

2 

1 

1 

3 

_ 4 

CUTOUT 

1 1 
12 

FOR  NUMBERING  OF  NODES 
ALONG  A  CROSS  SECTIONAL 
CIRCLE  NSCHM= 1 


NOTE:  Node  numbering  must  originate  at  the  upper  left-hand  corner  of  the 
finite-difference  grid. 


5.  PANEL  DIMENSIONS  CARD  (3F10.0) 
Columns  1-10 

11-20 


21-30 


XDIM  -  Length  of  a  panel  parallel  to 
a  generator. 

YDIM  -  Length  of  a  panel  along  a  cross 
sectional  circle  (not  the  pro¬ 
jection)  . 

R  -  Panel  radius  of  curvature.  (Any 
number  will  suffice  for  a  plate. 


6.  BOUNDARY  CONDITION  IDENTIFIER  CARD  (815) 

1 


Columns  1-5 
6-10 
11-15 
16-20 
21-25 
26-30 
31-35 
36-40 


NCASE1 

NCASE2 

NCASE3 

NCASE4 

NCASE5 

NCASE6 

NCASE7 

NCASE8 


=  <  3 


h 


4 


4 


4 
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NCASE_  identifies  a  boundary  of  the  panel  as  indicated  in  the  accompanying 
figure.  An  integer  from  1  to  5  is  assigned  to  each  NCASE_  accordinq  as  the 
desired  support  condition  along  the  edge  identified  by  NCASE_  is  one  of  those 
shown  above. 


NCASE1 


NCASE3 


D i rec t i on  of  a 


7.  RIGID  BODY  CONDITION  CARD  (415) 


Columns  1~5 

1  RBI 

6-10 

IRB2 

11-15 

IRB3 

16-20 

1 RB4 

Tangential  displacement  parallel 
to  the  boundary  identified  by 
I RB  is  precluded  at  the  mid¬ 
point  of  the  boundary  for  an 
even  number  of  grid  points,  or 
at  the  grid  point  nearest  mid¬ 
point  for  an  odd  number  of  qrid 
points. 


0  No  modifications  to  boundary 
conditions  are  made. 


8.  INITIAL-LOAD  CARD  (4F1 5 .0,15) 


Columns  1-15  PZSTRT  -  Initial-load  normal  to  the  panel 

surface.  (Positive  away  from  the 
center  of  curvature.) 

16-30  PXSTRT  -  Initial-load  parallel  to  a  genera¬ 

tor  of  the  panel.  (Positive  in 
the  direction  of  increasing  grid 
point  numbers.) 

31-45  PYSTRT  -  Initial-load  tangent  to  a  cross- 

sectional  circle.  (Positive  in 
the  direction  of  increasing  qrid 
point  numbers. 
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1)6-60 


XDISPO 


Prescribed  initial  axial  end- 
displacement.  Negative  for  com¬ 
pression. 


61-65 


MDISP 


10  Prescribed  loads 
1  Prescribed  end-displacements 


12 

23 

3** 


PXSTRT 

AND  PYSTRT  FOR  NSCHM=0 


zjJLS _ 

i 

i 

- 1 — —  PXSTRT 

PXSTRT 

POSITIVE  DIRECTIONS  FOR 


AND  PYSTRT  FOR  NSCHM=I. 


9.  LOAD  IDENTIFICATION  CARD  (1)15) 


Columns  1-5 


2 


LCASE  = 


3 


It 


Concentrated  load  at  grid 
point  number  LNDDE. 

Line- load  along  row  number  LROW 
of  the  finite-difference  grid. 

Line-load  along  column  number 
LCOL  of  the  finite-difference 
grid. 

Uniform  distributed  load  over 
the  entire  surface. 


6-10  LNODE  -  Number  of  the  finite-difference 

grid-point  at  which  a  concentrated 
load  i s  appl i ed . 

11*15  LROW  -  Number  of  the  row  of  finite- 

difference  grid  points  along  which 
a  line-load  is  applied. 

1 6—20  LCOL  *  Number  of  the  column  of  finite- 

difference  grid  points  along  which 
a  line-load  is  applied. 


139 


10.  LOAD  INCREMENT  CARD  (3F10.0) 


Columns  1-10 


2 1  NCR  -  Load  increment  normal  to  the  panel 
surface. 


11-20 


21-30 


XINCR  -  Load  increment  parallel  to  a 
generator  of  the  panel. 

YINCR  -  Load  increment  tangent  to  a  cross- 
sectional  circle  of  the  panel. 


Positive  directions  for  load- increments  are  the  same  as  those  for  the  initial 
loads . 

If  I  EX  t  1  cards  11,  12,  and  13  are  not  required. 


INITIAL  IMPERFECTION  CARD  (5F10.0.I5)  Required  only  if  LGRNG  =  0  and  I  EX 


Columns  1-10 


11-20 


21-30 


31-40 


W0  -  Amplitude  of  initial  geometric 
imperfect  ion. 

C0NS1  -  Wave  number  of  the  imperfection 
associated  with  the  direction  of 
a  generator  of  a  panel  (ir/half 
length  of  panel) 


C0NS2 


O  Clamped  curved  edges  and  unsup¬ 
ported  straiqht  edges. 

b  Clamped  curved  edges  and  simply- 
supported  straight  edges,  b  is 
one  half  of  the  circumferential 
length  of  a  panel 


X0  -  x-coordinate  of  geometric  center 
of  panel  measured  from  lower  left 
hand  corner 


1*1-50 


Y0  -  y-coord inator  of  geometric  center 
of  the  panel  measured  from  the 
lower  left  hand  corner 


51-60 


1MPF0RM  =■ 


Initial  imperfection  has  the  form 
W0(x,y)  -  WI*(l+cos  vt/a)  ( 1  - (n/b) 2 ) 
(Clamped  at  curved  edges  -  simply- 
supported  along  the  straight  edges) 

]  Initial  imperfection  has  the  form 
jW0(x)  =  WI*(l+cos  tt F, / a ) 

I  (Clamped  along  the  curved  edges  - 
[free  at  the  straight  edoes) 


12.  IMPERFECTION-GRID  PARAMETER  CARD  (215)  Required  only  if  LGRNG  =  1  and  I  EX  =  1 


Columns  1-5 


6-10 


MX  =  Number  of  grid  points  along  a 

generator  for  which  discrete  values 
of  initial  imperfection  are  known. 
Spacing  of  these  points  is  required 
to  be  uniform. 

MV  =  Number  of  grid  points  along  a  cross 
sectional  circle  for  which  discrete 
values  of  initial  imperfection  are 
known.  Spacing  of  these  points  is 
required  to  be  uniform. 


MX  or  MY  can  not  exceed  10  for  the  present  dimensions  of  the  program. 


13-  IMPERFECTION  DATA  CARD  (10F8.0)  Required  only  if  LGRNG  =  1  and  IEX  =  1 

A  single  card  contains  the  known  discrete  value  of  initial  imperfection 
for  every  grid  point  in  a  row  of  the  imperfection-grid.  Thus,  MY  cards 
are  required:  one  for  each  row  of  the  imperfect i on-qr i d .  A  row  is 
considered  to  be  parallel  to  a  generator  of  the  panel.  Imperfection 
values  are  positive  away  from  the  center  of  curvature  of  the  panel. 

Imperfection-grid  size  cannot  exceed  10  grid  points  in  either  direction 
unless  internal  dimensions  are  modified. 

Imperfection-grid  is  parallel  to  finite-difference  grid  with  origin  as 
shown  in  the  figure. 


Th  i  rd 
Second 
First 


D i rect i on  of  a 


Generator 


|1*.  NODAL  COORDINATE  CARD  (2F10.0)  Required  only  if  NAUTfc  =  0 

Columns  1-10  x-coordinate  of  a  finite-difference  grid 

point 

11-20  y-coordinate  of  a  finite-difference  grid 

point 

These  cards  are  required  only  when  NAUT0  =  0;  i.e.,  whenever  a  variable 
finite-difference  grid  is  required. 


11*1 


15.  STIFFENER  CARD  (if  15) 


Columns  1-5  NSTFR 

6-10  NSTYP 

11-15  NSLOC 

16-20  ISSPC 


/  0  no  s  t i f feners 
I  1  stiffeners  included 

J  1  through  8  -  cross  section  *ype 
'  (see  below) 

JO  located  on  outside  of  panel  (+) 
1  1  located  on  inside  of  panel 

JO  spaced  on  all  qrid  lines 
|  1  variable  spacing 


16.  STIFFENER  CROSS  SECTION  INPUT  CARDS  only  if  NSTFR  =  1 


The  number  and  required  input  of  these  cards  depend  on  the  choice  of 
stiffener  cross  section  type 


NSTYP  =  1  (2  FI  0.0)  (2  FI  0.0) 

(1)  Columns  1-10 

11-20 

(2)  Columns  1-10 

11-20 

NSTYP  =  2  (if  FI  0 . 0)  (2F10.0) 

(1)  Columns  1-10 

11-20 

21-30 

31-if0 

(2)  Columns  1-10 

11-20 

NSTYP  =  3  (5F10.0)  (2F10.0) 

(1)  Columns  1-10 

11-20 

21-30 

31-ifO 

ifl-50 

(2)  Columns  1-10 

11-20 


(NSTYP)  above. 

rectangu lar  plate 

HGHT  =  h 
W0TH  =  b 

SEM0D  =  elastic  modulus 
SGM0D  =  shear  modulus 

symmet  r i c  I  beam 

HGHT  =  h  (total  -  tf) 
WDTH  =  b 
FLTHK  =  tf 
WTHK  =  t 

w 

SEM0D  =  elastic  modulus 
SGM0D  =  shear  modulus 

non-symmetr i c  I  beam 

HGHT  =  h  (total  -  t,) 
TWDTH  =  b] 

BWDTH  =  b2 
FLTHK  =  tf 
WTHK  =  tw 

SEM0D  =  elastic  modulus 
SGM0D  =  shear  modulus 


J  I 


1  if  2 


J 


1 


NSTYP  =  4  (4F10.0) 

(2F10.0) 

inverted  T 

1 

( 1 )  Columns  1-10 

HGHT  =  h  (total  -  t,/2) 

f 

11-20 

WDTH  =  b  1 

21-30 

FLTHK  -  tf 

1 

31-40 

WTHK  =  tw 

r— -* 

-J 

- i-*t 

(2)  Columns  1-10 

SEM0D  =  elastic  modulus 

✓ 

lr 

11-20 

SGM0D  =  shear  modulus 

NSTYP  =  5  ( 4  F 1 0 . 0) 

(2F10.0) 

JT  beam 

! 

(1)  Columns  1-10 

HGHT  =  h  (total  -  tf/2) 

?  — — 1 

_ 1 

11-20 

WDTH  =  b  ! 

1 

21-30 

FLTHK  =  tf 

■i. 

31-40 

WTHK  =  tw 

J _ . 

_ 

(2)  Columns  1-10 

SEM0D  -  elastic  modulus 

7^ 

± 

w- 

11-20 

SGM0D  -  shear  modulus 

y 

NSTYP  =  6  (4F10.0)  (2FI0.0) 

( 1 )  Columns  1-10 

11-20 

21-30 

31-40 

(2)  Columns  1-10 

11-20 


C  beam 

HGHT  =  h  (total  -  tf) 
WDTH  =  b 
FLTHK  =  tf 
WTHK  =  tw 

SEM0D  =  elastic  modulus 
SGM0D  =  shear  modulus 


L 


j  't» 


r- 


V-t-1 


NSTYP  =  7  (4F10.0)  (2F10.0) 

(1)  Columns  1-10 

11-20 

21-30 

31-40 

(2)  Columns  1-10 

11-20 


L  beam 

HGHT  =  h  (total  -  tf/2) 
WDTH  =  b 
FLTHK  =  tf 
WTHK  =  tw 

SEM0D  =  elastic  modulus 
SGM0D  =  shear  modulus 


NSTYP  =  8  (4F10.0)  (5F10.0)  (2F10.0)  user's  choice 


r 


(1)  Columns  1-10  SAREA 

11-20  CEI 

21-30  ETAI 

31-40  ALPHA 


cross  section  area 

1st  principal  moment  of  inertia, 

2nd  principal  moment  of  inertia,  I,,., 
angle  from  horizontal  tog-direction, 
in  radians 


(2)  Columns  1-10 
11-20 
21-30 
31-40 
41-50 


CED  =  distance  to  contact  point,  ^-direction 
ETAD  =  distance  to  contact  point ,  7- di rect ion 
CES  =  distance  to  shear  center,  ^.direction 
ETAS  =  distance  to  shear  center,  w-direction 
C0NSTJ  =  twistinq  constant,  J  ^ 
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(3)  Columns  1-10 
11-20 


SEMOD  -  elastic  modulus 
SGMOD  -  shear  modulus 


Refer  to  figure  on  the  next  page  for  an  example  of  these  quantities 
for  the  user's  choice  stiffener,  type  8. 


NSTYP  =  5  (4F10.0,  15)  (5FI0.0)  QUAS  I  -  I  S0TR0P  I  C  HAT  SECTION 
(1)  Columns  1-10  HGHT  =  h 

WDTH  =  b 
FLNGW  =  b 
STHK  =  t 


(2)  Columns 


11-20 
21-30 
31  -  A0 
41  -  45 
1  -  10 


1 


NLAYER  =  No.  of 
layers. 

Ell  =  elastic 
modulus  parallel  to 
the  fibers  of  a  lamina. 


UN  IFORW  SffCT/6/V  TrttCKMESS 


11-20  E22  =  elastic  modulus 

perpend i cu 1 ar  to  the 
fibers  of  a  lamina. 


21  -  30  G12  =  lamina  shear  modulus 

31  -  40  PNU12  =  Poisson  ratio  associated 

with  stress  along  fiber  direction. 

41  -  50  PNU21  =  Poisson  ratio  associated 

with  stress  perpendicular  to  fiber 
d i rect i on . 


This  stiffener  cross  section  type  accomodates  f iber-reinforced 
stiffeners  where  the  fiber  directions  are  alternately  zero  and  90 
degrees.  Midplane  symmetry  is  assumed. 


Note  that  the  stiffeners  when  located  on  the  inside  of  the  panel  are 
inverted  as  depicted  below. 


5T/rre*i*R  orJ  ~rHe-  srir  ffajck.  o*j  tof. 

outs  tps  or  pa  net.  xrJsipe  of  FoUcl.  . 


NUMBER  OF  STIFFENERS  (15)  needed  only  if  ISSPC  =  1,  NSTFR  =  1 

Columns  1 -5  NST  =  number  of  stiffeners  if  they  are 

not  on  each  grid  line 


VARIABLE  SPACING  STIFFENER  CARDS  (15)  only  if  ISSPC  =  1,  NSTFR  =  1 

Columns  1-5  NSROW  =  row  number  that  locates  stiffener 

(if  NSCHM  =  0) 

column  number  that  locates  stiffener 
(if  NSCHM  =  1) 

One  card  is  required  for  each  stiffener,  so  there  will  be  the  same 
number  of  cards  as  NST  above.  The  row  (or  column)  numbers  must  be  in 
increasing  order,  i.e.,  1,**,7,10,  etc.  not  1,7, **,10. 

NUMBER  OF  LAYERS  CARD  (15) 

Columns  1_5  KN  -  Number  of  layers  in  the  laminate. 

LAYER  LOCATION  CARDS  (F10.0) 

Columns  1-10  Distance  (inches)  between  a  layer  surface 

nearest  the  center  of  curvature  and  the 
laminate  reference  surface.  First  layer 
is  considered  to  be  the  one  nearest  the 
center  of  curvature.  Distances  measured 
toward  the  center  of  curvature  are  negative. 

A  card  for  the  most  remote  surface  of  the  last  layer  is  required.  Thus, 
there  is  required  KN+1  LAYER  LOCATION  CARDS. 


21.  MATERIAL  PROPERTIES  CARDS  (6F10.0)  One  card  is  required  for  each  layer. 


Columns  1-10 

Modulus  of  elasticity  parallel  to  the 
fibers  of  a  given  layer  (lb/in2) 

11-20 

Modulus  of  elasticity  perpend i cu lar  to 
the  fibers  of  a  given  layer  (lb/in2) 

21-30 

Poisson  ratio  associated  with  strain 
parallel  to  the  fiber  axis  due  to  a  stress 
perpendicular  to  the  fiber  axis 

31-40 

Poisson  ratio  associated  with  strain  perpen¬ 
dicular  to  the  fiber  axis  due  to  a  stress 
parallel  to  the  fiber  axis. 

41-50 

Shearing  modulus  of  elasticity  (lb/in2) 

51-60 

Angle  between  the  fiber  direction  and  a 
generator  of  the  panel.  This  angle  is 
positive  whenever  the  fiber  direction  is 
situated  in  a  clockwise  orientation  relative 
to  the  positive  x-direction. 
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EXAMPLE  THAT  DEFINES  INPUT  QUANTITIES 
FOR  USERS  CHOICE  STIFFENER 


C  —  CSsTTZOiD 
S  —  SHtSAR  C&/STHR 

D  ~  CONTACT  Point 

BETW EEn  PA MEL 
ANO  STiFFBNBR. 
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